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ABSTRACT 



An impluse response h mn (q) of each transmission path is 
estimated from N received signals r m (m=l, . . . , M) and a 
known signal (for a number of users equal to N, n=l, . . . , 
N). MxN matrix H (q) having h mn (q) as an element and a 
QxQ matrix H having H(q) as an element are determined 
(where Q represents a number of multipaths of each trans- 
mitted wave and q=0, . . . , Q-l). A soft decision value b' n (k) 
is determined from decoded X 2 [b n (k)], and this is used to 
generate an interference component matrix B'(k) to generate 
an interference replica HB'(k). The interference replica 
HB'(k) is subtracted from a received matrix y(k) to deter- 
mine y'(k)- y(k) and H are used to determine an adaptive 
filter coefficient w n (k) to be applied to an n-th user in order 
to eliminate residual interference components in y'(k) 
according to the minimum mean square error criteria. y(k) is 
passed through w n (k) to provide a log-likelihood ratio as a 
received signal from the user n from which interferences are 
eliminated. 
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TURBO-RECEPTION METHOD AND 
TURBO-RECEIVER 

BACKGROUND OF THE INVENTION 

[0001] The invention relates to turbo -reception method 
and a turbo -receiver as may be used in a mobile communi- 
cation, for example, and which apply an iterative equaliza- 
tion utilizing a turbo -coding technique to waveform distor- 
tions which result from interferences. 

[0002] A task in the mobile station communication busi- 
ness is how to construct a system capable of acquiring a 
multitude of users on a limited frequency domain with a high 
quality. A multi-input multi-output (MIMO) system is 
known in the art as means for solving such a task. The 
architecture of this system is shown in FIG. 30A where a 
plurality of transmitters SI to SN transmit symbols c 1 (i) to 
c N (i) at the same time and on the same frequency, and the 
transmitted signals are received by an MIMO receiver 
equipped with a plurality of antennas #1 to #M. The received 
signals are processed by the receiver, which estimates trans- 
mitted symbols c 1 (i) to c N (i) from the respective transmitters 
SI to SN and delivers them separately as c-^i) to c N (i) to 
output terminals Out 1 to Out N. 

[0003] Up to the present time, a study of a specific 
implementation of an MIMO receiver in an MIMO system 
is not yet satisfactorily warranted. If one attempts to con- 
struct an MIMO receiver in an MIMO system on the basis 
of MLSE (maximum likelihood estimation) criteria, denot- 
ing the number of transmitters by N and the number of 
multi-paths through which a wave transmitted from each 
transmitter reaches the MIMO receiver by Q, the quantity of 
calculation required for the MIMO receiver will be on the 
order of 2 (Q " 1)N , and will increase even more voluminously 
with an increase in the number of transmitters N and the 
number of multi-paths Q. When information from a single 
user is transmitted as parallel signals, which are then 
received, a separation of individual parallel signals from 
each other requires a quantity of calculation, which increases 
exponentially with number of multi-paths. Accordingly, the 
present invention proposes herein an improved calculation 
efficiency turbo -reception method for a plurality of channel 
signals. To start with, an existing turbo-receiver for a single 
user (single transmitter) or a single channel transmitted 
signal, which illustrates the need for the present invention, 
will be described. 

[0004] Turbo-receiver for Single User 

[0005] An exemplary arrangement for a transmitter and a 
receiver is illustrated in FIG. 31. In a transmitter 10, 
information series c(i) is encoded in encoder 11, and an 
encoded output is interleaved (or rearranged) by an inter- 
leaver 12 before it is input to a modulator 13 where it 
modulates a carrier signal, the resulting modulated output 
being transmitted. The transmitted signal is received by a 
receiver 20 through a transmission path (each channel of 
multipath). In the receiver 20, a soft input soft output (SISO: 
single-input single output) equalizer 21 performs an equal- 
ization of delayed waves. At the input to the equalizer 21, the 
received signal is generally converted into a baseband, and 
the received baseband signal is sampled with a frequency 
which is equal to or greater than the frequency of symbol 
signals of information series in the transmitted signal to be 
converted into a digital signal, which is then input to the 
equalizer 21. 



[0006] For the single user, this corresponds to N=l in FIG. 
30A, and a received output form each reception antenna #m 
in (m=l, . . . , M) can be represented as follows: 



^(^)=2 q= o Q -^ m (^)^(^)+v m (A:) (1) 

[0007] where m represents an antenna index, h a channel 
value (a transmission path impulse response: a transmission 
path characteristic), b(k-q) a transmitted symbol from a user 
(transmitter 1), and v m (k) an internal thermal noise of the 
receiver 20. All outputs from the antennas #1 to #M are 
denoted by a matrix as indicated by an equation (2) to define 
an equation (3). 

K*M'i<*> 2 (*) ■ ■ ■ rUQT (2) 



^^Hiqybik-q^vik) 



(3) 



[0008] where 



v(*Hvi(*)v 2 (*) . - . v M (/:)] T (4) 

HiqyVKiq) . . . h M (q)Y (5) 

[0009] It is to be noted that [] T represents an inverted 
matrix. In cosideration of the number of channels Q of the 
mutipath, the following matrixes and matrix are defined: 



#Hr>e-i»e-2) 

[0010] where 



■(*)F 



(6) 
(7) 



// = 



(0) 



fiir>- 



(2-D 



0 



Hi 



(0) 



0 



(£-1) 



(8) 



b{k-q)=\b{k+Q-l)b{k+Q-2) . . . b(k-Q+l)y (9) 
n(£)=[v T (£+01)v T (£+0-2) . . . v T (A)] T (10) 

[0011] r(k) as defined above is input to the SISO equalizer 
21, which is a linear equalizer, deriving a log -likelihood 
ratio A 1 (LLR) of a probability that each encoded bit {b(i)} 
is equal to +1 to a probability that it is -1 as an equalization 
output. 



A:[W]=log 



Pr[b(k) = +l \y(k)] 
Pr[b(k) = -l \y(k)] 



(11) 



(12) 



[0012] where X 1 [b(k)] represents an extrinsic information 
fed to a succeeding decoder 24 and ^ 2 P [b(k)] a priori 
information applied to the equalizer 21. The log -likelihood 
ratio A 1 [b(k)] is fed to a sub tractor 22 where the a priori 
information >w 2 [b(k)] is subtracted therefrom. The result is 
then fed through a deinterleaver 23 to an SISO channel 
decoder 24, which calculates a log -likelihood ratio A 2 as 
follows: 



A 2 [*(i)]=log 



Pr[b(J) = +1 | A ![£>(/)], / = 0, ... ,B-1] 
Pr[b(i) = -1 | X l [£(/)], i = 0, ... ,5-1] 



(13) 



[0013] B: frame length 
^[fc(0]+V[K0] 



(14) 
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[0014] where X 2 [b(i) ] represents an extrinsic information 
which is applied as ^ 2 p [b(k)] to the equalizer 21 during the 
iteration, while ^i[b(k)] is applied as a priori information 
X 2 p [b(i)] to the decoder 24. In a subtracter 25, X 2 [b(i)] is 
subtracted from A 2 [b(i)], and the result is fed through an 
interleaver 26 to the equalizer 21 and the subtracter 22. In 
this manner, the equalization and the decoding are iterated to 
achieve an improvement of an error rate. 

[0015] To describe the prestage equalizer 21 in detail, the 
calculation of a linear filter response applied to a received 
matrix y(k) will be described. Using the a priori information 
X 2 p [b(k)] for the equalizer 21, a soft decision symbol esti- 
mate 



2/(£)=tan A[X 2 p [fc(£)]/2] 



(15) 



[0016] is calculated. Using the estimate and a channel 
matrix H, an interference component or a replica Hb'(k) of 
the interference component is reproduced and subtracted 
from the received signal. Thus, 



=H-{b{k)-b'{k))+n{k) 

[0017] where, 

b%k)=[b'(k+Q-l) . . . 0 . . . b'{k-Q+l)Y 



(16) 
(17) 

(18) 



[0018] Because the replica Hb'(k) of the interference 
component cannot be always a correct replica, the interfer- 
ence component cannot be completely eliminated by the 
equation (16). So a linear filter coefficient w(k) which 
eliminates any residue of the interference component is 
determined according to the MMSE (minimum mean square 
error) technique indicated below. 

w(A)=arg min \\^{k)yXk)-b(k)f (19) 

[0019] where H represents a conjugate transposition and | [ 
a norm. w(k) which mimimizes the equation (19) is deter- 
mined. 

[0020] Deriving w(k) in this manner is described in Daryl 
Reynolds and Xiandong Wang, "Low Complexity Turbo- 
Equalization for Diversity Channels" (http:/ee. tamu.edu/ 
Reynolds/). A major achievement of this technique lies in a 
significant reduction in the quantity of calculation. The 
quantity of calculation according to a conventional MLSE 
turbo has been proportional to the order of 2 Q_1 while a 
suppression to the order of Q 3 is enabled by this technique. 
It will be seen that w H (k)-y'(k) represents an output from the 
equalizer 21, and is used to calculate >. 1 [b(k)], which is then 
fed through the deinterleaver 23 to the decoder 24 to be used 
in the decoding calculation. 

[0021] For purpose of equalization in the equalizer 21, it 
is necessary to estimate the channel value (transmission path 
impulse response) h appearing in the equation (1). This 
estimation is hereafter referred to as a channel estimation. 
The channel estimation takes place by using a received 
signal of a known training series such as a unique word 
which is transmitted as a leader of one frame and a stored 
training series. A poor accuracy of the channel estimation 
prevents an equalization in the equalizer 21 from occurring 
in a proper manner. The accuracy of the channel estimation 
can be enhanced by increasing the proportion which the 
training series occupies in one frame, but this degrades the 
transmission efficiency of the intended data. Accordingly, it 



is desirable that the accuracy of the channel estimation could 
be improved while reducint the proportion of the training 
series in one frame. 

[0022] This is not limited to a receiver for multiple chan- 
nel transmitted signals inclusive of MIMO, but the same is 
true in the channel estimation of a receiver such as RAKE 
receiver or a receiver using an adaptive array antenna where 
the certainty of a decoded result is improved by an iterative 
decoding process. 

[0023] The described turbo -receiver has following restric- 
tions: 

[0024] It is an accommodation for a single user 
(single transmitter) or only one series transmitted 
signal. 

[0025] A channel value (matrix H) is necessary in 
reproducing an interference component, and this 
must be estimated in actual implementions. An esti- 
mation error results in a degradation in the effect of 
an iterative equalization. 

[0026] It is an object of the invention to provide compen- 
sations for these two restrictions by providing a turbo - 
reception method and a receiver therefor which allow the 
receiver mentioned above to be extended to a receiver for a 
plurality of transmitted series signals such as for multiple 
users or parallel transmissions from a single user. 

[0027] It is another object of the invention to provide a 
reception method and a receiver therefor in which a channel 
value of a received signal is estimated from the received 
signal and a known signal serving as a reference signal, the 
received signal is processed using the estimated channel 
value, and the processed signal is decoded so that the 
processing using the estimated channel value and the decod- 
ing are iterated upon the same received signal and which 
allow the channel estimation to be achieved with good 
accuracy using a relatively short known signal. 

SUMMARY OF THE INVENTION 

[0028] According to a first aspect of the present invention, 
there is provided a turbo-reception method of receiving N 
series (where N is an integer equal to or greater than 2) 
transmitted signals. The method comprises calculating a 
channel value h mn (q) (n=l, . . . , N) from M received signals 
r m (m=l, . . . , M) and N series known signals, determining 
a soft decision transmitted symbol b' n (k) on the basis of N 
series a priori information ^ 2 [b n (k)] which are obtained by 
the decoding, and using the channel value h mn (q) and the 
soft decision transmitted symbol b' n (k) to calculate an inter- 
ference component H-B'(k) formed by an intersymbol inter- 
ference produced by an n-th transmitted signal itself and 
transmitted signals other than the n-th transmitted signal as 
follows: 



H = 



Hi 



(0) ■•■ H(Q-D 



0 



(0) 



0 



(e-D 
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-continued 



B'(/r)=[fc' T (*+G-l) . . • . . . Z/ T (£-e+l)] T 

&t£+?H&'i(*+?)&' 2 (fc+<?) ■ ■ ■ b' N (k+q)Y 
q=Q-l . . . -Q+l for q*Q 
b'{k)=[b\{k) . . . 0 . . . b r ^(k)Yq=0 

[0029] where b'(k) has a zero element at n-th position, Q 
represents the number of multipaths of each transmitted 
signal wave, q=0, . . . , Q-l and [] T represents a transposed 
matrix, subtracting the intersymbol interference HB'(k) 
from the received matrix y(k) to obtain a difference matrix 
y'(k) where 

y(JfcHr T (*K2-l>- T (*H2-2) ■ ■ ■ r^(k)Y 

[0030] determining an adaptive filter coefficient W n (k) 
applied to the received signal for the n-th transmitted signal 
in order to eliminate residual interference component in the 
difference matrix y'(k) using the channel matrix H or refer- 
ence signal, filtering the difference matrix y'(k) according to 
the adaptive filter coefficient W n (k) to obtain a log-likeli- 
hood ratio for the n-th series as the received signal for the 
n-th series transmitted signal from which the interference 
has been eliminated, and decoding by using the the log- 
likelihood ratios for N series. 

[0031] According to a second aspect of the present inven- 
tion, in the arrangement according to the first aspect fo the 
invention, when q=0, 

b'{k)=[b\{k) . . . -Ab'M) . . . b'^k)T 

[0032] where an element f(b' n (k)) is located at n-th posi- 
tion, and f( ) represents a function which satisfies f(0)=0 and 
having as a variable b' n (k) which satisfies d{f(b' n (k))}/ 
d{b' n (k)}^0. 

[0033] According to a third aspect of the present inven- 
tion, the equalization takes place in a plurality of stages, and 
the number of series in the equalization output is sequen- 
tially reduced in the later stages. 

[0034] According to a fourth aspect of the present inven- 
tion, in a turbo-reception method in which a channel value 
of a received signal is estimated from the received signal and 
a known signal serving as a reference signal, the received 
signal is processed using the estimated channel value, per- 
forming a decoding processing upon the processed signal 
and the processing which uses the estimated channel value 
and the decoding processing are iterated upon the same 
received signal, the certainty that the decoded hard decision 
information signal has is determined on the basis of the 
value of an associated soft decision information signal, and 
a hard decision information signal having a certainty which 
is equal to or greater than a given value is also used as a 
reference signal in the channel estimation of the next itera- 
tion. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0035] FIG. 1 is a diagram showing a functional arrange- 
ment of a system including an embodiment of a turbo- 
receiver according to a first aspect of the present invention; 



[0036] FIG. 2 is a schematic view showing a specific 
example of a functional arrangement of a multiple output 
equalizer 31 shown in FIG. 1; 

[0037] FIG. 3 is a flow chart of an embodiment of a 
turbo-reception method according to the first aspect of the 
present invention; 

[0038] FIG. 4A is a diagram of an exemplary frame; 

[0039] FIG. 4B is an illustration of processings which are 
performed during each iteration, presented in order to illus- 
trate an iterative channel estimation method according to a 
fourth aspect of the present invention. 

[0040] FIG. 5 is a diagram of an exemplary functional 
arrangement for deriving a hard decision symbol which is 
likely to be certain; 

[0041] FIG. 6 is a flow chart of an exemplary processing 
procedure for channel estimation according to the present 
invention; 

[0042] FIG. 7A is a diagram of an exemplary functional 
arrangement of part of the equalizer 31 according to the 
second aspect of the invention which reflects an error 
correcting decoding result of a signal being detected; 

[0043] FIG. 7B is a diagram showing an exemplary pro- 
cessing procedure therefor; 

[0044] FIG. 8 is a diagram of an example of a receiver 
which uses a turbo -equalizer iteratively; 

[0045] FIG. 9 is a diagram of an exemplary receiver 
which performs RAKE reception-turbo-decoding itera- 
tively; 

[0046] FIG. 10 is a diagram of an exemplary receiver 
which performs an adaptive array antenna reception- turbo - 
decoding iteratively; 

[0047] FIG. 11A is a schematic view of a turbo -equalizer; 

[0048] FIG. 11B is a schematic view of a turbo -decoder; 

[0049] FIG. 12 is a schematic view of a receiver which 
iterates a processing of a received signal using an estimated 
channel and a decoding processing of the processed signal; 

[0050] FIG. 13 is a flow chart of an exemplary processing 
procedure of a reception method which iterates a processing 
of a received signal using an estimated channel and a 
decoding processing of the processed signal; 

[0051] FIG. 14A is a diagram of an exemplary frame 
arrangement; 

[0052] FIG. 14B is a schematic illustration of an iterative 
processing of an estimation of a channel H and a noise 
covariance matrix U when a received signal contains noises 
other than white Gaussian noise; 

[0053] FIG. 15 is a schematic view of an exemplary 
functional arrangement of part of an equalizer which is used 
in the estimation of the noise covariance matrix U; 

[0054] FIG. 16 is a flow chart of an exemplary processing 
procedure which iterates the estimation of a channel value 
used in the estimation of the noise covariance matrix U and 
the decoding processing; 
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[0055] FIG. 17 is a diagram showing the principle of a 
turbo -receiver according to a third aspect of the present 
invention; 

[0056] FIG. 18 is a diagram of an exemplary functional 
arrangement of a turbo-receiver according to the third aspect 
of the present invention; 

[0057] FIG. 19 is a diagram showing a specific arrange- 
ment of a functional arrangement of a multi-user (prestage) 
equalizer 71 shown in FIG. 18; 

[0058] FIG. 20 is a flow chart of an exemplary processing 
procedure of a turbo -reception method according to the third 
aspect of the present invention; 

[0059] FIG. 21 is a diagram showing another exemplary 
functional arrangement of a multi-stage equalization section 
according to the third aspect of the present invention; 

[0060] FIG. 22 is a diagram of an example of a system 
arrangement in which an embodiment according to the first 
aspect of the present invention (2) is applied; 

[0061] FIG. 23 graphically shows an error rate character- 
istic of a turbo -receiver according to the first aspect of the 
present invention (1), with assumption that the channel has 
been perfectly estimated, E b representing a bit power and N Q 
a noise power; 

[0062] FIG. 24 graphically shows an error rate character- 
istic when the iterative channel estimation has been per- 
formed while changing a threshold value (Th); 

[0063] FIG. 25 graphically shows an error rate character- 
istic of a turbo -receiver according to the fourth aspect of the 
present invention, in particular, employing the iterative 
channel estimation; 

[0064] FIG. 26 graphically shows an error rate character- 
istic of the turbo-receiver which employs an estimation of a 
noise co variance matrix U; 

[0065] FIG. 27 graphically shows an error rate character- 
istic of the turbo -receiver shown in FIG. 1; 

[0066] FIG. 28A graphically shows an error rate charac- 
teristic as plotted against E b /N Q of an embodiment according 
to the second aspect of the present invention which reflects 
an error correcting decoding result of a signal being 
detected; 

[0067] FIG. 28B graphically shows an error rate charac- 
teristic plotted against a of the embodiment according to the 
second aspect of the present invention; 

[0068] FIG. 29 graphically shows a result of a simulation 
of an error rate characteristic of a turbo-receiver according 
to the third aspect of the present invention; 

[0069] FIG. 30A is a diagram illustrating the concept of a 
MIMO system; 

[0070] FIG. 30B is a diagram of an arrangement in which 
received signals from a pair of antennas are input to a 
turbo -receiver as four received signal series; and 

[0071] FIG. 31 is a diagram of a functional arrangement 
of a conventional turbo -transmitter and receiver for a single 
user. 



EMBODIMENTS OF THE INVENTION 

[0072] First Aspect of the Invention (1) 

[0073] FIG. 1 shows an exemplary arrangement of an 
MIMO system to which the present invention is applied. 

[0074] In each of N transmitters SI . . . SN, information 
series c 1 (i) to c N (i) are encoded in encoders 11-1, . . . , 11-N, 
and the encoded outputs are fed through in terle avers 12-1, . 
. . , 12 -N to modulators 13-1, . . . , 13-N as modulation 
signals, thus modulating a carrier signal in accordance with 
these modulation signals to transmit signals b 1 (k) to b N (k). 
In this manner, transmitted signals b 2 (k) . . . b N (k) from the 
transmitters SI, . . . , SN form N series transmitted signals. 

[0075] A received signal r(k) which is received by a 
multiple output receiver through transmission paths (chan- 
nels) is input to a multiple output equalizer 31. A signal 
received by the receiver is converted into a baseband signal, 
which is then sampled at one-half the symbol period ,for 
example, to be converted into a digital signal, which is then 
input to the equalizer 31. It is assumed that there are one or 
more digital signals, the number of which is represented by 
an integer M. For example, received signals from M anten- 
nas are formed into received signals in the form of M digital 
signals. 

[0076] The equalizer 31 delivers N log-likelihood ratios 
Aifb^k)], . . . A-^b^k)]. In subtracters 22-1, . . . , 22-N, a 
priori information X 1 [b 1 (k)], . . . ^ 1 [b N (k)] are subtracted 
from A 1 [b 1 (k)], . . . A x [b N (k)], and results are fed through 
deinterle avers 23-1, . . . , 23 -N to be input to soft input soft 
output (SISO) decoders (channel decoders) 24-1, . . . , 24-N 
to be decoded, whereby the decoders 24-1, . . . , 24-N deliver 
decoded information series c\(i) . . . c' N (i) together with 
A 2 [b 1 (i)], . - . , A 2 [b N (i)]. In subtractors 25-1, . . . , 25-N, 
X 1 [b 1 (i)] ? . . . >^ 1 [b N (i)] are subtracted from A 2 [b 1 (i)], . . . , 
A 2 [b N (i)], respectively, and results are fed through interleav- 
ers 26-1, . . . , 26-N to be fed as ^ 2 [ D i(X)] ? • • • XJb^k)] to 
the multiple output equalizer 31 and the subtractors 22-1, . 
. . , 22-N, respectively. 

[0077] A received signal r m (k) (m=l, . . . , M) from 
multiple users (a plurality of transmitters) is a sum of 
received signals from the plurality of users when it is input 
to the equalizer 31 as indicated below. 



1VV ^q=0 ^n=l 



te>K(k-q)+vjjc) (20) 

[0078] where q=0, . . . , Q-l, and Q represents the number 
of multipaths for each transmitted wave. Defining a matrix 
y(k) by a similar procedure as used for the single user, we 
have 



=H-B(k)+n(k) 

[0079] where 



(21) 
(22) 



[0080] 



r 



M 



(0) 



0 



H (0) 



0 



H {Q-l) 



(23) 



where 
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-continued 



^(9) = 



^11(<?) ■■■ h\N(q) 



h 



Ml{q) - - - hMN(q) 



(24) 



B(k)=[b\k+Q-l) . . . 6 T (£) . . . b T (k-Q+l)J 

b(k+q)=[b 1 (k+q)b 2 (k+q) . . . b N (k+q)Y q=Q-l, Q-2, . 
■ , -Q+l 



(25) 



(26) 



[0081] In an interference elimination step, it is assumed 
that it is desirable to obtain a signal from an n-th user 
(transmitter). In this example, a soft decision symbol esti- 
mate for signals from all users (transmitters) and a channel 
matrix (transmission path impulse response matrix) H are 
used to produce a synthesis of interferences by user signals 
other than from n-th user and interferences created by the 
n-th user signal itself, or an interference replica H B (k) is 
reproduced. Then, the interference replica is subtracted from 
y(k) to produce a difference matrix y'(k). 



(27) 
(28) 

(29) 



y '(k)=y(k)-H-B \k) 
=H-(B{k)-B'{k))+n{k) 

[0082] where 

B(k)=[b T (k+Q-l) . . . b fT (k) . . . b'^(k-Q+l)T 

b'{k+q)=[b\{k+q)b'?(k+q) . . . b' N (k+q)] T : q=Q-l, . . . 

, -Q+l, q*0 (30) 

ft'(*Hb'i(ifc) - ■ ■ 0 . . . b'^(k)Y: q=0 (31) 

[0083] b'(k) has a zero element at n-th position. It is to be 
understood that b f n (k) represents a soft decision transmitted 
symbol estimate which is obtained by calculation of b' n (k)= 
tan h[X 2 [b n (k)]/2]. The matrix B'(k) represents a replica 
matrix of the interference symbol. 

[0084] A filter coefficient w n (k) for the n-th user, which is 
used to eliminate the residue of the interference component, 
namely, a residual interference based on the imperfectness of 
the interference component replica H-B'(k) and an interfer- 
ence component created by the n-th signal itself, is deter- 
mined according to MMSE (minimum mean square error) 
criteria as one which minimizes the following equation (32): 



w n (£)=arg min|K H (*)y(*)-^(*)|f 



(32) 



[0085] Subsequent operation remains the same as for the 
single user. Specifically, w n (k) which is obtained in this 
manner is used to calculate w n H (k)-y'(k), and a result of 
calculation is fed through deinterleaver 23-n to be input as 
^i[b n (i)] to the decoder 24-n where a decoding calculation 
is made. 

[0086] The described method of applying a filter (linear 
equalization) processing upon the received signal r m is 
repeated for all of users 1 to N. As a consequence, the 
number of outputs from the equalizer 31 will be equal to N, 
and all these outputs are decoded by respective decoders 
24-1, . . . , 24-N. What has been mentioned above is an 
extension of a turbo-receiver for single user to a receiver for 
multiple users (MIMO). 

[0087] From the foregoing description, it will be seen that 
an exemplary functional arrangement of the multiple output 
equalizer 31 will be as shown in FIG. 2. Specifically, M 
received signals r m (k) are supplied to a received matrix 
generator 311 where a received matrix y(k) is generated and 



is fed to equalizers 312-1 to 312 -N for each user. The 
channel matrix H which is calculated by a channel estimator 
28 is also supplied to the equalizers 312-1 to 312-N. A priori 
information ^J^nOO] from each channel decoder 24-n is 
supplied to a soft decision symbol estimator 313 where a soft 
decision transmitted symbol estimate b' n (k)=tan h[X 2 [b n (k)]/ 
2] is calculated. All of equalizers 312-1 to 312-N have 
identical functional arrangements and process in an identical 
manner, and accordingly, a typical one of them (312-1) will 
be described. 

[0088] In addition, soft decision transmitted symbol esti- 
mates b' n (k) to b' N (k) are supplied to an interference replica 
matrix generator 314-1 where an interference replica matrix 
B' 1 (k) is generated according to the equations (29) to (31), 
then the matrix B' 1 (k) is subject to a filter processing 
according to the channel matrix H in a filter processor 315-1, 
and a resulting interference replica component H-B' 1 (k) is 
subtracted from the received matrix y 2 (k) in a difference 
calculator 316-1 to produce a difference matrix y , 1 (k). 

[0089] At last the channel matrix H or a reference signal 
which will be described later is input to a filter coefficient 
estimator 317-1 to determine the filter coefficient w 2 (k) 
which is used to eliminate the residue of the interference 
component. In the example shown, the channel matrix H and 
a covariance a 2 of a noise component and soft decision 
transmitted symbols b' 1 (k) to b' N (k) from the soft decision 
symbol generator 313-1 are input to the filter coefficient 
estimator 317-1 to determine the filter coefficient w 1 (k) 
which minimizes the equation (32) according to the mini- 
mum mean square error criteria. A specific example of 
determining the filter coefficient w 1 (k) will be described 
later. The difference matrix y'i(k) is processed with the filter 
coefficient w 1 (k) in an adaptive filter processor 318-1, and 
A 1 [b 1 (k)] is delivered as an equalized output for the received 
signal which corresponds to the transmitted signal from user 
1. 

[0090] A processing procedure for the multiple input mul- 
tiple output turbo -reception method according to the 
described embodiment of the present invention is shown in 
FIG. 3. At step SI, a channel value h mn (q) and a covariance 
a 2 of a noise component are calculated from a received 
signal r(k) and each training signal b n (k). At step S2, the 
channel matrix H is calculated from the channel value 
h mn (q). At step S3, the soft decision transmitted symbol 
estimate b' n (k)=tan h[X 2 [b n (k)]/2] is calculated on the basis 
of the a priori information X 2 [b n (k)] which is obtained 
during a previous iteration of the turbo-reception processing. 

[0091] At step S4, the received signal y(k) is generated 
from the received signal r(k). At step S5, the interference 
replica matrix B' n (k) is generated according to the equations 
(29) to (31) using the respective soft decision transmitted 
symbol estimates b' n (k) . At step S6, the interference com- 
ponent replica HB' n (k) for the received signal from the n-th 
transmitter is calculated. At step S7, the interference com- 
ponent replica H-B' n (k) is subtracted from the received 
matrix y(k) to provide the difference matrix y' n (k). At step 
S8, using the channel matrix H, the soft decision transmitted 
symbols b' 1 (k) to b' N (k) and the covariance a 2 of the noise 
component, the filter coefficient w n (k) which is used to 
eliminate the residue of interference remaining in the 
received signal from the n-th transmitter is determined 
according to the minimum mean square error criteria so as 
to minimize the equation (32). 
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[0092] At step 9, a filter processing according to the filter 
coefficient w n (k) is applied to the difference matrix y' n (k) to 
obtain the log- likelihood ratio A 1 [b n (k)]. At step S10., the 
a priori information X 2 [b n (k)] is subtracted from A 1 [b n (k)] 
and the result is deinterleaved and decoded to deliver the 
log-likelihood ratio A 2 [b n (k)]. The steps S4 to S10 are 
performed either simultaneously or sequentially for n=l to 
N. Subsequently, at step Sll, an examination is made to see 
whether or not the number of decoding operations, namely, 
the number of turbo-reception processings has reached a 
given number. If the given number is not reached, the 
operation goes to step S12 where the extrinsic information 
X 2 [b n (k)] is subtracted from the log-likelihood ratio 
A 2 [b n (k)], and its result is interleaved to determine the a 
priori information ^ 2 [b n (k)], thus returning to step S3. If it 
is found at step Sll that the number of decoding operations 
has reached the given number, a resulting decoding output is 
delivered at step S13. 

[0093] The channel estimator 28 will now be described. 
Each received signal r m (k) can be represented as follows: 

r^^^-^^h^iqyb^k-q^v^k) (33) 

[0094] The channel estimator 28 determines the channel 
value (transmission path impulse response) h mn (q) appear- 
ing in the equation (33) and the mean power (so 2 ) of the 
noise v m (k). Normally, a unique word (training signal) 
which is known to the receiver is inserted at the beginning 
of each frame to be transmitted on the transmitting side, as 
shown in FIG. 4 A, and the receiver utilizes the unique word 
(known signal) as a training series to estimate the channel 
value h mn (q) using RLS (recursive least square) technique. 
Each of the channel decoder 24-1, . . . , 24-N delivers +1 if 
a log-likelihood ratio A 2 [b 1 (i)], . . . A 2 [b N (i)] is positive, and 
-1 if the latter is negative as a decoded code signal (or 
transmitted encoded symbol hard decision value) 6 1 (i), . . . 
, 6 N (i), which are fed iteratively through interleavers 27-1, 
. . . , 27-N to the channel estimator 28. The received signal 
r(k) is input to the channel estimator 28 together with a 
unique word which is fed from a unique word memory 29 to 
serve as a reference signal. On the basis of these input 
signals, the channel estimator 28 estimate each h mn (q) 
according to the equation (33) and a 2 value according to the 
least square technique. This estimation can be made in a 
similar manner as the estimation of an impulse response 
when adaptively equalizing the received signal with an 
adaptive filter by estimating an impulse response of a 
transmission path. 

[0095] The use of the training series in this manner is 
conventional, but in order to enhance the net transmission 
rate, it is necessary to reduce the proportion that the unique 
word occupies in one frame, but this increases the error of 
the channel estimation. If there is such an error, it results in 
degrading the iterative equalization response mentioned 
above. According to the present invention, an iterative 
estimation of the channel value in the following manner is 
proposed. 

[0096] The concept according to the present invention is 
illustrated in FIG. 4B. The purpose of the concept is to 
estimate the channel value iteratively in each stage of the 
iterated equalization processing of the same received signal 
or during the each iteration of the turbo -reception process- 
ing. During the first run, for information symbol series 
which follows the unique word, the channel value is esti- 



mated using only the unique word as a reference signal, and 
the estimated channel value is used to equalize the received 
signal and to estimate the transmitted symbol. However, 
before the second iteration of the equalization processing, 
the channel estimation is made using the unique word as the 
reference signal, and the symbol estimate (hard decision 
value) which is obtained during the previous decoding 
processing is also used as the reference signal to perform a 
channel estimation of the entire frame. In this instance, 
rather than using every hard decision value, only that hard 
decision value or values which have been determined as 
likely to be certain is used as the reference signal. The hard 
decision is rendered by delivering +1 when the logic- 
likelihood ratio A 2 [b n (i)] from the decoder 24-n is positive, 
and delivering -1 if the later is negative. The greater the 
absolute magnitude of the log-likelihood ratio A 2 [b n (i)], it 
can be said the greater the likelihood that the hard decision 
value is likely to be certain. For example, the certainty that 
is rendered when the log-likelihood 5 is determined to be +1 
is greater than the certainty when the logic-likelihood 0.3 is 
determined to be +1. Accordingly, an interative channel 
estimation method which uses a threshold value to designate 
hard decision value b n (i) which is likely to be certain will be 
described below. 

[0097] Initially, using the log -likelihood ration A 2 [b n (i)] 
from the decoder 24-n, a soft decision symbol value b' n (i) is 
determined as follows: 

fr'n(0=tan A[A 2 [6 n (0I/2] 

[0098] This operation is made in order to normalize the 
logic-likelihood value to 1, thus preventing its absolute 
magnitude from exceeding 1. Next, a threshold between 0 
and 1 is provided, and when the absolute magnitude of a soft 
decision value b' n (i) is greater than the threshold value, a 
corresponding hard decision value £> n (i) is preserved to be 
used in the iterative channel estimation. For example, if the 
threshold value is chosen to be 0.9, those of hard decision 
values 6 n (i) which correspond to soft decision values b' n (i) 
having absolute magnitudes equal to or greater than 0.9 are 
selected. It is believed that the certainty of the selected hard 
decision value 6 n (i) is high in view of the fact that the 
threshold value is as high as 0.9, and accordingly, it is 
believed that the accuracy of the iterative channel estimation 
which is made using such hard decision values can be 
improved. However, it is also considered that because the 
number of selected symbols is reduced, the accuracy of the 
iterative channel estimation may be influenced and becomes 
degraded. Accordingly, it is necessary that an optimum 
threshold value be choosed between 0 and 0. To add, if the 
threshold value is to be 1, if follows that there is no hard 
decision value o n (i) which is selected, stating to the effect 
that no iterative channel estimation takes place. Accordingly, 
a threshold value on the order of 0.2 to 0.8 is chosen, as will 
be further described later. 

[0099] Accordingly, those of transmitted symbol estimates 
(hard decision values) 6 1 (i), . . . , D N (i) for the information 
symbol series during the first transmission which are deter- 
mined to be likely to be certain according to the threshold 
value are fed from the outputs of the interleavers 27-1, . . . 
, 27-N to a previous symbol memory 32 to be stored therein 
as a previous transmitted symbol estimate. During the sec- 
ond iteration of the equalization and decoding processing for 
the received signal r(k) (it being noted that the received 
signal r(k) is stored in a memory), the unique word is 
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initially used to make a channel estimation, and those of 
estimated transmitted symbol hard decision estimates D 1 (i), 
. . . , 6 N (i) which have been determined to be likely to be 
certain are read from the previous symbol memory 32 to be 
fed to the channel estimator 28 to make a channel estima- 
tion, namely, a channel estimation for the entire frame. A 
resulting estimate h mn (q) and a 2 are used to perform an 
equalization and decoding (the estimation of transmitted 
symbol) with respect to the received signal r(k). At this time, 
those symbol values among the estimated transmitted sym- 
bols which had been determined to be likely to certain 
according to the threshold are iused to update the content of 
the previous symbol memory 32. Subsequently, an estima- 
tion using the unique word and an estimation using those of 
the previously estimated transmitted signals which are deter- 
mined to be likely to be certain are used to make a channel 
estimation of the entire frame during a subsequent iteration 
of the equalization and decoding. The estimated channel is 
used to perform the equalization and decoding (the estima- 
tion of transmitted symbol) and to update the previous 
symbol memory 32. Alternatively, those of the transmitted 
symbol hard decision values 6 1 (i), . . . , o N (i) from the 
decoders which are determined to be likely to be certain 
according to the threshold value may be directly stored in the 
previous symbol memory 32 to update it, and when the 
symbol values stored in the previous symbol memory 32 are 
to be utilized, they will be fed through the interleavers 27-1, 
. . . , 27-8 to be input to the channel estimator 28. 

[0100] By the iterations which proceed in this manner, an 
error of the channel estimation can be reduced, the accuracy 
of the symbol estimation can be improved and the problem 
of a degraded response in the turbo-equalization due to an 
error of an channel estimation can be improved. 

[0101] When a channel estimation is made with respect to 
information symbol series using the symbol hard decision 
values which are likely to be certain in the manner men- 
tioned above, a function as indicated in FIG. 5 is added to 
each decoder 24-n. The log- likelihood ratio A 2 [b n (i)] is 
input to a soft decision value estimator 241 to calculate 
b' n (i)=tan h[A 2 [b n (i)]/2] thus estimating a transmitted sym- 
bol soft decision value b' n (i). The value b' n (i) is compared 
against a threshold value Th from a threshold pre setter 243, 
thus delivering 1 when b' n (i) is equal to or greater than Th 
and delivering 0 when it is less than Th. On the other hand, 
the log-likelihood ratio A 2 [b n (i)] is input to a hard decision 
unit 244, which delivers +1 if A 2 [b n (i)] is positive and 
delivers -1 if the latter is negative as a symbol hard decision 
value D n (i). The symbol hard decision value D 1 (i) is fed 
through a gate 245 which is opened when a corresponding 
symbol soft decision value is equal to or greater than the 
threshold value, and is passed through the interleaver 27 -n 
shown in FIG. 1 to be fed to the previous symbol memory 
32, thus updating the stored symbols. 

[0102] A procedure of the channel estimation which also 
utilizes a symbol hard decision value or values which are 
likely to be certain is shown in FIG. 6. Initially, at step SI, 
a channel estimation is made with a received signal r(k) and 
a unique word. At step S2, an examination is made to see if 
the decoding processing is for the first time, and if it is, at 
step S3, the estimated channel value h mn (q)is used to per- 
form the equalization and decoding processing or the opera- 
tions shown at steps S3 to S10 in FIG. 3. 



[0103] At step S4, a transmitted symbol hard decision is 
rendered with respect to the log-likelihood ratio A 2 [b n (i)] to 
determine a hard decision value 6 n (i). At step S5, b' n (i)=tan 
h[A 2 [b n (i)]/2] is calculated from the log -likelihood ratio 
A 2 [b n (i)] , thus estimating a transmitted symbol soft decision 
value b' n (i). At step S6, by examining whether the symbol 
soft decision value b' n (i) is equal to or greater than the 
threshold value Th or not, those of the symbol hard decision 
values 6 n (i) which are likely to be certain are determined. At 
step S7, the symbol hard decision value or values which are 
likely to be certain are used to update the content of the 
previous symbol memory 32. At step S8, an examination is 
made to see whether or not the number of decoding opera- 
tions has reached a given value, and if not, the operation 
returns to step SI, or more exactly, going through the step 
S12 shown in FIG. 3 and returning to the step SI shown in 
FIG. 3. 

[0104] If it is found at step S2 that the decoding processing 
is not for the first time, at step S9, a previous stored symbol, 
namely, a hard decision symbol which is likely to be certain 
is read from the previous symbol memory 32, and it is used 
together with information symbol series of the received 
signal r(k) to make a channel estimation, thus transferring to 
step S3. 

[0105] In the above description, even during a second and 
a subsequent iteration, the channel estimation is made from 
the beginning using the unique word as a reference signal, 
but it is also possible that only the hard decision symbol 
which is likely to be certain may be used as the reference 
signal during a second and a subsequent iteration. In this 
instance, as indicated in broken lines in FIG. 6, an exami- 
nation is made to see if this processing is for the first time 
at step SI', and if it is, the unique word is used as a reference 
signal together with the received signal to estimate the 
channel value. After storing the estimated channel value and 
parameter values which are used in the estimation in a 
memory at step S3', the operation transfers to the equaliza- 
tion and the decoding processing which take place at step 
S3. 

[0106] If it is found at step SI' that the processing is not 
for the first time, before the channel estimation is made, a 
previously stored channel estimate and various processing 
parameters are set up at S4', whereupon the operation 
transfers to step S9. 

[0107] It is to be noted that a solution of the equation (32) 
is given as follows: 

w n (A:)=(//G : (>t)// :H +a 2 /)- 1 -/? (34) 

[0108] where I represents a unit matrix, a 2 an internal 
noise power of a receiver (covariance of noise component), 
a 2 I a covariance matrix of noise component and G(k) 
corresponds a squared error of the channel estimation. 

G(k)=E[ (B(k)-B '(£)) • (B(k)-B '(£)) H ]=diag[D(£+e- 1 ), . 

. . ,D(k), . . - , D(k-Q+iy] (35) 

[0109] where E [] represents a mean, diag a diagonal 
matrix (having elements which are not along the diagonal 
being zero's). 

D(k+q)=dmg[l-b a 1 (k+q), . . . , l-b'\{k+b), . . . , 

l-b%(k+q)-] (36) 

q=Q-l, Q-2, . . . , -Q+l, q*0 
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[0110] and when q=0, 

D(k)=di ag [l-b'\(k), 1-6' 2 N (£)] (37) 

[0111] i Y appearing in the matrix D(k) represents an n-th 
element (it is assumed that an n-th user's transmitted signal 
is assumed to be a desired signal). 

(38) 

h = . 



[0112] Thus, h comprises all the elements in the (Q-1)N+ 
n-th column of H appearing in the equation (23). Input to the 
filter coefficient estimator 317-1 of the multiple output 
equalizer 31 as shown in FIG. 2 are the channel matrix H 
and the noise power a 2 which are estimated in the channel 
estimator 28 and the soft decision transmitted symbols b' n (k) 
to b' N (k) from the soft decision symbol generator 313-1 to 
calculate the residual interference elimination filter coeffi- 
cient w n (k) according to the equations (34) to (38). 

[0113] It will be seen that the equation (34) requires an 
inverse matrix operation, but the required calculation can be 
reduced by using the matrix inversion lemma for the inverse 
matrix. Specifically, by approximating each b' 2 appearing in 
the equations (36) and (37) by 1, it follows that 

D(k+q)=dmg[0, . . . , 0]=0(^0) (39) 
D(*)=diag[0, . . . , 1, . . . , 0] (40) 

[0114] Thus, D(k) have elements having a value of 1 for 
those elements which are located at n-row and n-column 
while all other elements are equal to 0. When the error 
matrix G(k) of the equation (35) which is determined by the 
equations (39) and (40) is substituted into the equation (34), 
we have 

w n (^)=(A-A H +a 2 /)- 1 -/j (41) 
[0115] where h is as defined by the equation (38). 

[0116] With this approximation, w n (k) does not depend on 
k, and accordingly, an inverse matrix calculation at every 
discrete time k can be dispensed with, thus reducing the 
quantity of calculation. 

[0117] The matrix inversion lemma for the inverse matrix 
is applied to the equation (41). The lemma states that 
assuming A and B are (M, M) square matrices, C a (M, N) 
matrix and D a (N, N) square matrix, when A is given by 
A _1 =B _1 +CD _1 C H , the inverse matrix of A is given as 
follows: 

A- ± =BBC{D+C™BBC)- 1 C^B (42) 

[0118] Applying the lemma to the inverse matrix operation 
appearing in the equation (41), we have 

hikyhik^+c^^B^+CD- 1 ^ 

hikyhik^CD-^tO^B^thik^C 
/= J D" 1 ,A(yt) H =C H 

[0119] Using these equations to calculate the equation 
(42), the inverse matrix operation appearing in the equation 
(41) can be calculated. While the equation (42) contains an 
inverse matrix operation (D+C H BBC) _1 , this inverse matrix 
becomes scalar and can be simply calculated. 



[0120] Thus, in this instance, it is reduced to the following 
form: 

w n (Jt)=l/(o 2 +/7 H -/z> (41-1) 

[0121] l/( ) on the right side of this equation is scalar or 
assumes a constant value, which may be chosen to be 1. 
Thus, we can put w n (k)=h, indicating that w (k) can be 
determined by only h. As indicated in broken lines in FIG. 
2, only h which is represented by the equation (38) in the 
channel matrix H from the channel estimator 28 may be 
input to the filter coefficient estimator 317-1. 

[0122] Without being limited to the use of the matrix 
inversion lemma for the inverse matrix for the approxima- 
tion by the equations (39) and (40), the approximation also 
allows the quantity of calculation for the equation (34) to be 
reduced. In particular, when this approximation is used, and 
the matrix inversion lemma for the inverse matrix is used, 
the quantity of the calculation can be further reduced. In 
addition, if the covariance matrix of the noise component is 
represented by a 2 I, an approximation: w n (k)=h may be used 
as indicated by the equation (41-1), whereupon it is inde- 
pendent from the covariance matrix, allowing a further 
simplification of the calculation. 

[0123] Second Aspect of the Invention (reflecting an error 
correction) 

[0124] In the equalization processing where H-B'(k) is 
subtracted from the received matrix y(k) represented by the 
equation (27), an error correction decoding result is reflected 
in a transmitted signal soft decision value for a signal other 
than the signal b n (k) being detected, but an error correction 
decoding result which relates to the signal b n (k) being 
detected is not. For this reason, it is preferred to employ the 
following processing. 

[0125] b'(k) appearing in the equation (29) or the equation 
(31) is changed as follows; 



b>(K)=[b\{k)b'^k) . . . b'^ikyfib'Jkftb^k) . . . 



(43) 



[0126] where f(b' n (k)) is an arbitrary function having 
b' n (k) as an input. 

[0127] When such change is made, it becomes possible to 
reflect an error correction decoding result with respect to the 
signal b n (k) being detected. Thus, rather than using b' n (k)=0, 
by adding a suitable value depending on b' n (k), it is possible 
to emphasize a signal being detected which is buried in a 
noise or an interfering signal, thus allowing b n (k) to be 
properly detected. 

[0128] Because the sign of b' n (k) relates to a result of a 
hard decision rendered upon a symbol which corresponds to 
b' n (k), and in view of the fact that the greater the absolute 
magnitude of b' n (k), the higher the reliability of the hard 
decision symbol which corresponds to b' n (k), it is necessary 
that f(b' n (k)) satisfies the following requirements: 

[0129] The function f has a value of 0 for b' n (k)=0 or when 
the reliability of the hard decision symbol is equal to 0. Or, 



/(0)=0 



(44) 



[0130] In addition, the greater the value of b' n (k), the 
greater the value of the function f. Or, 



d{f(b' n (k))}/d{b' n (k)}^0 



(45) 
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[0131] Examples of such f(b' n (k)) include the following: 

/(fc' n (*))=axb' n (*) (46) 
/(£' n (£))=axb' n (A;) 2 (47) 

[0132] For example, when the equation (46) is used and a 
is chosen to be a constant, the equation (43) can be imple- 
mented in a simple manner. Here, a must satisfy the 
requirement: 0<a<0.6. If a is greater than 0.6, BER (error 
rate) will be degraded, preventing a correct decoding result 
from being obtained. It is also contemplated to make a to be 
variable in accordance with the reliability of a decoding 
result. For example, a. may be chosen for each iteration of 
the decoding processing. In this instance, the reliability of 
the decoding result will be normally improved with the 
number of iterations of the decoding processings and 
accordingly, an increased value may be chosen for a depend- 
ing on the number of iterations of the decoding processing. 
Alternatively, the reliability of an entire decoded frame may 
be determined upon each iteration of the decoding process- 
ing, and the value of a may be chosen in accordance with the 
reliability thus determined. To determine the reliability of 
the decoded frame, a decoding result may be compared 
against a decoding result which is obtained during an 
immediately preceding iteration, and a number of hard 
decision symbols which changed from the previous decod- 
ing operation may be counted, for example. Thus, where 
there is an increased number of changed decision symbols, 
the reliability may be determined to low while when the 
number of changed hard decision symbols is small, the 
reliability may be determined to be high. 

[0133] As b' n (k) is changed in this manner, the equation 
(35) which is used when determining the filter coefficient 
w n (k) for MMSE (minimum mean square error) filter is 
desirably changed as follows: 

G(^=£[(^(/r)-5'(^))-(5(A:)-^'(^)) H ]=diag[Z)(^+C?-l), . 
. . , D(k), D(k-Q+iy\ 

[0134] Using the equations (29) and (31), it follows that 
assuming 

B'(k)=[bXk+Q-l)b'(k+Q-2) . . . b'{k) . . . b'(k-Q+l)Y 
b\k)=[b'Jk)b^(k) . . . -f(b' n (k)) . . . b'^(k)F 
b'{k+q)=[b> ± {k+q)b'^k+q) . . . -f(b ! Jk+q)) . . . b'^(k+ 

q)Y:q*Q, q=Q-% ■ ■ ■ , -G+i 

[0135] D(k) has elements located at n-row and n-column, 
which are represented as follows: 

£[(fcn(^)+/(*'n(^)))-(^)+/(i'n(^)))*] 

[0136] where ( ) represents a complex conjugate. For 
BPSK modulation, this expression turns into the following 
expression: 

£[fo n ^) 2 + 2fe n W(fo n (A:))+/(fo n (A:)) 2 ]=£[6 n 2 (A;)] + 2£- 
[b n (k)f(b n (k)) + E[f(b n (ky) 

[0137] The first term has a mean value of 1. When b n (k) 
is approximated by b'(k), the equation (37) turns into the 
following form: 

D{k)=di^[l-b a 1 {k)l-b a - z {k) . . . l-b'\^(k) 1+2E 
V(b' n (k)b' n tk)) + E\f{h' n (kf)l-b'\ +1 {k) . . . l-b'\(k)\ (48) 

[0138] For example, when the equation (46) is chosen for 
f(b' n (k)), D(k) turns into the following form: 

Z)(A)=diag[l-^' 2 1 (A:)l-fo' 2 2 (A) . . . l-Z/ 2 n _ 1 (£)l+(2a+ 
a^b r2 n (k)l-b ,2 n+l(k) . . . l-b'\(ky\ (49) 

[0139] An exemplary functional arrangement which esti- 
mates an adaptive filter coefficient w n (k) when reflecting an 



error correction decoding result into a signal being detected 
is shown in FIG. 7A where a signal being detected is chosen 
as a transmitted signal b 1 (k) from a first transmitter. A soft 
decision transmitted symbol b' 1 (k) is input to a function 
calculator 331-1 to calculate a function f^b'-^k)). Soft deci- 
sion transmitted signal b f 1 (k) to b' n (k) from N decoders and 
f(b , 1 (k)) are input to an error matrix generator 332-1 to 
calculate and generate an error matrix G(k) according to the 
equations (35), (36) and (48). The error matrix G(k), an 
estimated channel matrix H and noise power a 2 are input to 
a filter coefficient generator 333-1 where a calculation of the 
equation (34) is made to estimate an adaptive filter coeffi- 
cient w n (k). In this instance, f(b' n (k)) is also input to an 
interference replica matrix generator 314-1, thus generating 
an interference replica matrix B'(k) represented by the 
equation (29) according to the equations (30) and (43). The 
filter coefficient W n (k) is used in filtering a difference matrix 
y'(k) in an adaptive filter 318-1, thus yielding a log-likeli- 
hood ratio A 1 [b 1 (k)]. It will be noted that in the filter 
coefficient estimator 317-1 shown in FIG. 2, the function 
calculator 331-1 shown in FIG. 7A is omitted, and only the 
soft decision transmitted symbols b , 1 (k) to b' N (k) are input 
to the error matrix generator 332-1 in order to calculate the 
equation (34). 

[0140] In the flow chart of FIG. 3, the interference replica 
matrix B'(k) is generated at step S4, and after processings at 
steps S5 to S7, the filter coefficient w n (k) is determined at 
step S8. If a calculation of the equation (34) is made during 
the processing at step S8, the soft decision transmitted 
symbols b' 1 (k) to b' N (k) are used to calculate the equations 
(35) to (37) to generate an error matrix G(k) at step S8-2, and 
the error matrix G(k), the estimation channel matrix H and 
the noise power are used to calculate the equation (34) to 
determine an adaptive filter coefficient w n (k) at step S8-3, as 
shown in FIG. 7B. 

[0141] When it is desired to reflect an error correction 
decoding result in a signal being detected as mentioned 
previously, a soft decision transmitted symbol b' n (k) of a 
signal which is to be detected may be calculated at step S8-1 
before entering the step S4, and this may be used at step S4 
where the equation (31) is replaced by the equation (43), or 
in other words, the equations (29), (30) and (43) may be used 
to generate an interference replica matrix B'(k), and at step 
S8-2, the equation (37) may be replaced by the equation 
(48). In the event f(b' n (k)) is chosen to be equal to a.b' n (k) 
or ctb' n (k) 2 and when a is chosen to be variable, a may be 
determined in accordance with the number of processing 
operations or the reliability of the entire decoded frame at 
step S8-1-1, and l+(2a+a 2 )b' n (k) 2 may be calculated and 
used as f(b' n (k)) at step S8-1-2. 

[0142] The technique of reflecting an error correcting 
result into a signal being detected is also applicable to a 
single user turbo -receiver which has been described initially 
in connection with the prior art. In the technique of reflecting 
an error correcting result into a signal being detected, the 
approximation represented by the equations (39) and (40) 
can be applied, and in this instance, only a matrix h shown 
by the equation (38) which is supplied from the channel 
estimator 28 may be input to the filter coefficient generator 
333-1, as indicated in broken lines in FIG. 7A. 

[0143] In the above description, the adaptive filter coeffi- 
cient w n (k) is determined according to the equation (34) or 
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by using the channel matrix H, but the use for the channel 
matrix H can be dispensed with. Specifically, during the 
initial decoding processing (turbo -reception processing), the 
error matrix G appearing in the equation (34) becomes a unit 
matrix. Accordingly, the difference matrix y'(k) and the 
training signal either alone or in combination with a hard 
decision transmitted symbol 6 n (k), preferably D n (k) having a 
high reliability in the sense mentioned above are input to the 
filter coefficient generator 333-1 to calculate the adaptive 
filter coefficient w n (k) in a sequential manner by application 
of RLS (recursive least square) technique. Because the error 
matrix G depends on a discrete time k, during a second and 
a subsequent iteration of the decoding operation, it becomes 
necessary to update the adaptive filter coefficient w n (k) from 
symbol to symbol, and as mentioned previously, it is pre- 
ferred to determine the adaptive filter coefficient w n (k) by 
using the channel matrix H. 

[0144] Fourth Aspect of the Invention (channel estima- 
tion) 

[0145] Using not only known information such as a unique 
word in the iterative channel estimation, but also using a 
hard decision value of information symbol, in particular one 
which is likely to be certain as a reference signal is appli- 
cable not only in the described multiple input multiple 
output turbo -reception method, but also generally to a turbo- 
reception method in which a channel (transmission path) of 
a received signal is estimated from the received signal and 
the known signal, the estimated channel value is used to 
process the received signal and to decode it, and the decoded 
signal is used in iterating the processing according to the 
estimated channel value and the decoding processing upon 
the same received signal. 

[0146] FIG. 8 shows an example in which a hard decision 
value of information symbol is used in a channel estimation 
turbo -equalizer 41. The turbo-equalizer 41 determines a 
linear equalization filter coefficient in accordance with an 
estimated channel value. The received signal is processed by 
such a linear equalization filter, the processed signal is 
decoded, and the decoded signal is utilized in the iterative 
processing of the same received signal. A received signal 
r(k) is input to the turbo-equalizer 41 and is also fed to a 
channel estimator 42 where a channel value (transmission 
path characteristic) is estimated on the basis of the received 
signal r(k) and a unique word from a memory 29. The 
received signal r(k) is subject to an equalization processing 
in the turbo -equalizer 41 in accordance with the estimated 
channel value, and is then subject to a decoding processing, 
whereby decoded data c'(i) and a soft decision value b'(i) are 
delivered. The soft decision value b'(i) is input to a symbol 
designator 43. If its soft decision value b'(i) has an absolute 
magnitude which is equal to or greater than a threshold, the 
corresponding hard decision value 6(i) is stored in a previ- 
ous symbol memory 32 as one which is likely to be certain 
(as having a high reliability) to update it. During a subse- 
quent iterative reception processing (equalization process- 
ing) of the same received signal r(k), the channel estimation 
which takes place in the channel estimator 42 takes place by 
using not only the unique word but also the hard decision 
value B(i) of information symbol which is stored in the 
previous symbol memory 32. 

[0147] The turbo -equalizer 41 may comprise the receiver 
shown in FIG. 1, for example, from which the iterative 



channel estimator 28, the unique word memory 29 and the 
previous symbol memory 32 are removed. It may comprise 
a receiver shown in FIG. 31. Again, a solution of the 
equation (19) will be in the following form according to the 
Wiener solution: 

[0148] where H is as defined by the equation (8), and 

h=[H(Q-l), . .. , H(0)Y 

[0149] where H( ) is as defined by the equation (5), 
a 2 =E[||v|| 2 ] (variance of noises), and 

A(k)=diag[l-b' 2 (k+Q-l), . . . , 1, . . . , l-fc' 2 (k-Q+l)] 

[0150] In this manner, also in the receiver shown in FIG. 
31, the channel matrix H( ) is estimated, and the channel 
matrix H( ) is used to determine the equalization filter 
coefficient w(k), the received signal is filtered according to 
the filter coefficient w(k), and the processed output is subject 
a decoding processing. Accordingly, by using a hard deci- 
sion information symbol of a high reliability in the channel 
estimation during the iterative recessive processing, there 
can be obtained a more correct channel estimation. 

[0151] FIG. 9 shows an example of a turbo -receiver in 
which the iterative channel estimation method is applied in 
the iterative reception in which RAKE synthesis takes place. 
A received signal r(k) is fed to a RAKE synthesis processor 
45 and a channel estimator 42. During an initial reception, 
a channel value is estimated on the basis of the received 
signal r(k) and a unique word in the channel estimator 42, 
and a compensation for a phase rotation which each symbols 
has experienced on a transmission path as well as RAKE 
synthesis take place in accordance with the estimated chan- 
nel value in the RAKE synthesis processor 45, or a time 
diversity processing takes place to be delivered to a turbo - 
decoder 46. The turbo -decoder 46 delivers decoded data c'(i) 
and a soft decision value b r (i). The soft decision value b'(i) 
is input to a symbol designator 43, and as in described 
examples, a hard decision value 6(i) of an information 
symbol which is likely to be certain is stored in a previous 
symbol memory 32 by updating it. During a second and a 
subsequent iterative reception processing of RAKE recep- 
tion turbo-decoding, not only the unique word but also the 
hard decision value of information symbol which is obtained 
during the previous iteration are used in the channel esti- 
mator 42 for purpose of channel estimation, thereby making 
the channel estimation more accurate to improve the quality. 

[0152] FIG. 10 shows an example of a turbo-receiver in 
which the described iterative channel estimation method is 
applied in the iterative reception using an adaptive array 
antenna. A received signal r(k) is received by an adaptive 
array antenna receiver 47, and is then branched into a 
channel estimator 42 where a channel estimation takes place 
in accordance with the received signal in combination with 
a unique word. The estimated channel value is used to 
determine the weights to be applied to each antenna element 
or corresponding reception paths in an array weight deter- 
mining unit so that the principal beam of the antenna 
directivity response of the adaptive array antenna receiver 
47 is directed toward the oncoming direction of an intended 
wave while null is directed toward the oncoming direction of 
an interfering wave, and such weights are applied to appli- 
cable locations. A received output from the adaptive array 
antenna receiver 47 is fed to a turbo-decoder 46 to be 
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decoded, whereby decoded data c'(i) and a soft decision 
value b'(i) are delivered. The soft decision value b'(i) is input 
to a symbol designator 43, and a hard decision value which 
is likely to be certain is stored in a previous symbol memory 
32 by updating it. During a second and a subsequent iterative 
reception processing of the adaptive array antenna receiver 
47-turbo-decoder 46, the hard decision value of information 
symbol which is obtained during the previous iteration is 
used in the channel estimation in the channel estimator 42 
together with the unique word. In this manner, the channel 
estimation takes place more correctly, resulting in a more 
accurate control of the antenna directivity response and 
improving the quality. 

[0153] The turbo-equalizer 41 shown in FIG. 8 is sche- 
matically shown in FIG. 11A, and as shown, it comprises a 
series connection of a soft input soft output (SISO) equalizer 
41a and an SISO decoder 41b, and the operation iterates 
between the equalizer 41a and the decoder 41b. The turbo- 
decoder 46 shown in FIGS. 9 and 10 is schematically 
shown in FIG. 11B, and as shown, it comprises a series 
connection of a SISO decoder 46a and an SISO decoder 466, 
and an iterative decoding operation takes place between the 
decoders 46a and 46b. The turbo-decoder 46 shown in 
FIGS. 9 and 10 may comprise a single SISO decoder. 

[0154] The examples shown in FIGS. 8 to 10 are collec- 
tively illustrated in FIG. 12. Thus, a received signal is 
initially processed in an iterative receiver (turbo -receiver) 49 
in accordance with a channel value which is estimated by a 
channel estimator 42, the processed signal is decoded, and a 
result of the decoding operation is provided in the form of 
decoded data (symbol) c'(i) and its soft decision value b'(i). 
The soft decision value b'(i) is compared against a threshold 
value in a symbol designator 43 to determine whether or not 
a corresponding decoded data c'(i) (symbol hard decision 
value) is likely to be certain. If it is determined to be likely 
to be certain, the hard decision value is stored in a previous 
symbol memory 32 by updating it, and during a second and 
a subsequent iteration of the processing using the estimated 
channel value and the decoding processing, the symbol hard 
decision value obtained during the previous iteration is used 
in addition to known information such as unique word in the 
channel estimation which takes place in the channel estima- 
tor 42 to provide a more correct channel estimation. 

[0155] FIG. 13 shows an exemplary processing procedure 
of an iterative turbo-reception method which also uses a 
symbol hard decision value. At step SI, a channel value is 
estimated on the basis of a received signal and a known 
signal. At step 2, an examination is made to see if this is the 
first one of the iterative processing, and if it is, the channel 
value which is estimated at step SI is used to process the 
received signal at step S3, and then the decoding processing 
takes place to determine a symbol hard decision value and 
a soft decision value. At step S4, a symbol hard decision 
value which corresponds to the symbol soft decision value 
and which is likely to be certain is extracted, and at step S5, 
the extracted symbol hard decision value is used to update 
a previous symbol hard decision value which is stored in the 
memory 32. At a step S6, an examination is made to see if 
the number of decoding processings has reached a given 
number, and if not, the operation returns to step SI. When 
it is found at step S2 that the iterative processing is not for 
the first time, a previous symbol hard decision value is read 
out from the memory 32 at step S7, and it is used together 



with an information symbol of a received signal in perform- 
ing a channel estimation, whereupon the operation transfers 
to step S3. 

[0156] Again, as mentioned previously in connection with 
steps SI' to S4' with reference to FIG. 6, the known signal 
may not be used during a second and a subsequent iteration. 

[0157] In the example shown in FIG. 10, RAKE synthesis 
processor 45 may be inserted between the adaptive array 
antenna receiver 47 and the turbo -decoder 46, as indicated 
by broken lines. In this instance, a channel estimation which 
is necessary for the compensation of a phase rotation of each 
symbol and the RAKE synthesis in the RAKE synthesis 
processor 45 may be served by the channel estimator 42, or 
may take place separately. 

[0158] Noise other than White Gaussian Noise 

[0159] In the embodiment of the turbo-reception method 
(according to the first aspect of the invention), the embodi- 
ment according to the second aspect of the invention which 
takes an error correction into consideration and the embodi- 
ment of the turbo-reception method characterized in its 
channel estimation method (according to the fourth aspect of 
the invention), the processing took place on an assumption 
that the noise is white Gaussian noise. Specifically, v m (k) 
appearing on the right side of the equation (29) indicating a 
received signal r m (k) from each antenna is assumed to be 
white Gaussian noise. What is meant by white Gaussian 
noise is a signal which follows the Gaussian distribution and 
have statistical features expressed as follows: 

( cr 2 for q = 0 



[0160] where E[] represents an expected value and a 2 a 
variance. White Gaussian noise may be exemplified by 
thermal noise which is generated in an antenna element. 
What is influenced by the assumption of white Gaussian 
noise is a 2 I portion appearing in the equation (34) which 
determines the filter coefficient w n (k) or the equation (50) 
which determines the filter coefficient w n (k). For example, 
w n (k) appearing in the equation (34) is calculated through 
the process indicated below. 

w n (/t)=(//G(A)i/ H +£[n(/t)-w H (yt)])- 1 /j=(//G(/t)// H +a 2 /)- 
ih 

[0161] where v m (k) is calculated as E[n(k>n H (k)]=a 2 I by 
the assumption of the white Gaussian noise having a vari- 
ance 6. The channel matrix H which is estimated by the 
iterative channel estimator 28 (FIG. 1) or 42 (FIG. 12), a 2 
and the error matrix G(k) which is calculated on the basis of 
the a priori log -likelihood value are substituted into the 
equation (34) to calculate the filter coefficient w n (k). 

[0162] When the noise v m (k) is not white Gaussian noise, 
E[n(k)-n H (k)]=a 2 I does not apply. Accordingly, in order to 
calculate the filter coefficient w n (k), it is necessary to 
estimate an expected value (covariance) matrix 
E[n(k)n H (k)] for the noise component by a separated 
method. Such method will now be described. A covariance 
matrix for the noise component will be abbreviated as 
U=E[n(k) n H (k)]. y(k)=H B(k)+n(k) in the equation (22) is 
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modified into n(k)=y(k)-HB(k) and is substituted into the 
covariance matrix U, as indicated below. 

U=E[n(kyn M (k)]=E[(y(k)-H-B(k))-(y(k)-H-B(k)) M ] 

[0163] If we can assume that a matrix y(k) is available 
from a received signal, an estimate H of a channel matrix H 
from the channel estimate and B(k) is available from a 
reference signal, it is possible to estimate the matrix U 
according to the time average method as follows: 

U=^o TT (y{k)-H-B (it)) ■ 0 (Jfc) -H-B{kyf (5 1) 

[0164] where Tr represents the number of reference sym- 
bols. 

[0165] During an iterative channel estimation which takes 
place in the iterative channel estimator 28 or 42, the cova- 
riance matrix U is estimated using the channel matrix H 
together with the equation (51). A procedure therefor is 
illustrated in FIG. 14. FIG. 14A shows a unique word and 
information symbol series which occurs in one frame of a 
received signal, and FIG. 14B shows the initial processing 
as well as a subsequent processing. During the initial pro- 
cessing, only the unique word is used as a reference signal, 
initially estimating the channel matrix H. Subsequently, U is 
estimated using the unique word and the channel matrix 
estimate H according to the equation (51). Using the esti- 
mate U and H, the filter coefficient w n (k) is calculated as 
follows: 

wJk^HGifyH^+lf)-^ (52) 

[0166] and the filter coefficient w n (k) is used to apply a 
first equalization upon the received signal, thus estimating 
transmitted information symbol. 

[0167] During a second iteration, the unique word as well 
as one of information signals (*) estimated during the initial 
equalization which is determined to be likely to be certain 
according to the threshold value are both used as reference 
signals to estimate H again with the same procedure as used 
during the initial processing, thus estimating U again. As this 
operation is repeated, the channel matrix estimate H 
becomes more accurate with the iteration, and the estimate 
of U becomes more accurate, thus improving the accuracy of 
the filter coefficient w n (k) to improve the response of the 
equalizer. 

[0168] In this manner, a turbo-reception when a noise 
other than white Gaussian noise is contained in a received 
signal is made possible. 

[0169] A functional arrangement in which a linear equal- 
ization is performed by estimating a covariance matrix U of 
a noise contained in a received signal is shown in FIG. 15 
as an example in which a log-likelihood ratio A 1 [b 1 (k)] is 
obtained as an equalization output from a multiple output 
equalizer 31 shown in FIG. 2 for a received signal corre- 
sponding to a transmitted signal from a first transmitter. 
Parts corresponding to those shown in FIG. 2 are designated 
by like reference numerals in FIG. 15 as used in FIG. 2. 

[0170] A unique word from a unique word memory 29 or 
a previous symbol hard decision from a previous symbol 
memory 32 which is likely to be certain is input to a 
reference matrix generator 319, which then generates a 
reference matrix B(k) according to the equations (25) and 
(26). The reference matrix B(k), an estimation channel 
matrix H from a channel estimator 28, and a received matrix 
y(k) from a received matrix generator 311 are supplied to a 



covariance matrix estimator 321, which then calculates the 
equation (51) to obtain an estimated matrix U for a covari- 
ance matrix U. 

[0171] Soft decision transmitted symbols b 1 (k), . . . , b n '(k) 
from a soft decision symbol generator 313-1 are input to an 
error matrix generator 322-1, where an error matrix G 1 (k) 
corresponding to the square error of the channel estimation 
is generated according to the equations (35), (36) and (37). 
The error matrix G 1 (k) 7 the estimated covariance matrix U 
and the estimatied channel matrix H are supplied to a filter 
estimator 323-1, where the equation (52) is calculated to 
estimate a filter coefficient w 2 (k). The filter coefficient w 1 (k) 
and the difference matrix y'(k) from a difference calculator 
316-1 are fed to an adaptive filter 318-1 where a filter 
processing w 1 (k) H y'(k) is applied to y'(k), and its result is 
delivered as a log-likelihood ratio A 1 [b 1 (k)]. 

[0172] When reflecting an error correction recording 
result into a signal being detected, a function calculator 
331-1 as shown in FIG. 7A is provided, as indicated by 
broken lines in FIG. 15, to calculate f(b' n (k)), and an 
interference replica matrix generator 314-1 uses the equation 
(43) rather than the equation (31), and a error matrix 
generator 322-1 uses the equation (48) rather than the 
equation (37). 

[0173] The procedure shown in FIG. 14B is shown as a 
flow chart in FIG. 16. Specifically, at step SI, a received 
signal r(k) and a known signal (for example, a unique word) 
are used to estimate a channel matrix H, and at step S2, an 
examination is made to see if this processing is for the first 
time, and if it is, the known signal, the estimated channel 
matrix H' and the received signal r(k) are used to calculate 
the equation (51) to determine an estimated covariance 
matrix U at step S3. 

[0174] At step S4, the estimated channel matrix H, the 
estimated covariance matrix U and an error matrix G(k) 
which comprises symbol soft decision values are used to 
calculate the equation (52) to estimate a filter coefficient 
w n (k). 

[0175] At step S5, the estimated channel matrix H and the 
filter coefficient w n (k) are used to equalize the received 
signal or to calculate the equation (27) to determine 
w n H (k)-y'(k) to obtain a log-likelihood ratio A-^b^k)], sub- 
sequently performing a decoding process to estimate a hard 
decision value and a soft decision value of a transmitted 
symbol. 

[0176] Purpose of step S6 is to determine a symbol hard 
decision value which corresponds to a symbol soft decision 
value which exceeds a threshold value and which is likely to 
be certain (or having a high reliability). This symbol hard 
decision value is used to update a symbol hard decision 
value which is stored in a previous symbol memory 32. 
Subsequently, an examination is made at step S8 to see if the 
number of the decoding processings has reached a given 
value, and if not, the operation returns to step SI. However, 
if a given number is reached, the processing upon this 
received frame is completed. 

[0177] If it is found at step S2 that the iterative processing 
is not for the first time, namely, for a second or a subsequent 
iteration, a symbol hard decision value is read from the 
previous symbol memory 32 at step S9, and is used together 
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with information symbol in the received signal to estimate 
the channel matrix H, subsequently transferring to step S3. 

[0178] Again, by changing the steps SI and 52 in the 
similar manner as steps SI' to S4' shown in broken lines in 
FIG. 6, it is possible to avoid the use of the known signal 
during a second and a subsequent iteration. When it is 
desired to reflect an error correction recording result into a 
signal being detected, a function f(b' n (k)) may be calculated 
at step S10, as shown in broken lines in FIG. 16, and it may 
be used to obtain the error matrix G(k). In either instance, 
the hard decision transmit symbol may not be used in the 
estimation of the covariance matrix U. The capability to 
estimate a covariance matrix U of a noise contained in a 
received signal which is other than white Gaussian noise is 
applicable to a variety of useful purposes, as will be 
described below. 

[0179] (1) A reception method for a multiple series trans- 
mitted signal containing an unknown interfering signal is 
cited. As shown in FIG. 30A, it is assumed that an inter- 
fering signal i(k) which is unknown to a turbo -receiver such 
as signals from other cells or zones in a mobile communi- 
cation, for example, to be received by the turbo -receiver in 
addition to N series transmitted signal as may be transmitted 
from N user transmitters. In this instance, the equation (20) 
can be written into the following form: 

[0180] In this model, putting i(k)+v m (k)=v' m (k), we have 

rJfc)=2 q . 0 ^/U?)-fcn(H)+v' m (^) (20) 1 ' 

[0181] Treating v' m (k) as a noise signal other than white 
Gaussian noise, H and U are estimated in a manner men- 
tioned previously, and w n (k) is estimated, and a turbo- 
reception can be made by iterating an equalization of a 
received signal and an estimation of transmitted symbols. 

[0182] (2) In a communication system which employs a 
transmission/reception separation filter, when an oversam- 
pling of a received signal is made at a higher rate than a 
symbol period, there occurs a correlation between noise 
components which are contained in received signals which 
are sampled at individual times, and this prevents the noises 
in the received signals from being treated as white Gaussian 
noise. In other words, the equation (20) does not apply. 
Accordingly, an assumption represented as 

E[n(k)-n H (k)]=a 2 I 

[0183] does not hold. A processing upon a received signal 
which is separated by 

[0184] the transmission/reception separation filter may 
utilize equation (51) to determine a covariance matrix U, 
thereby allowing the received signal to be properly pro- 
cessed. 

[0185] (3) In the described turbo-reception method, the 
arrangement is such that every multipath component from 
each transmitter (user) on Q paths are synthesized. However, 
in the event there exists a prolonged delay wave on channels 
(for example, assume that paths include one symbol delay, 
two symbol delay and three symbol delay path and there 
exist separately a thirty symbol delay path: in this instance, 
the thirty symbol delay path component is treat as a pro- 
longed delay wave), it is possible to prevent the prolonged 
delay wave from being synthesized, but to treat it as an 
unknown interference which can be eliminated by an adap- 



tive filter. When the prolonged delay wave component is 
treated as the interfering signal i(k) in the example according 
to the first aspect of the invention (1), it may be eliminated. 

[0186] In the processing of a received signal containing a 
noise other than white Gaussian noise, the estimation of the 
covariance matrix U is applicable to a single user turbo- 
reception method by allowing it to serve in place of a 2 ! in 
the equation (50). In a similar manner, it may be used in a 
RAKE synthesis processing reception illustrated in FIG. 9 
or a turbo -reception using an adaptive array antenna recep- 
tion shown in FIG. 10, irrespective of a single user or 
multi-user application, or more generally, in a channel 
estimation in a channel estimator 42 and an estimation of the 
covariance matrix U during the iterative decoding operation 
shown in FIG. 12. For RAKE reception, only the channel 
estimation may be used. 

[0187] Third Aspect of the Invention (multistage equal- 
ization) 

[0188] In the forgoing description, received signals r 1? . . 
. , r m are equalized in a multiple output equalizer 31 to 
determine log -likelihood ratios A 2 [b(k)], . . . , A N [b(k)], but 
in a modification (2) of the first aspect of the present 
invention, there are provided a plurality of equalizer stages 
in cascade connection in a manner such that the number of 
outputs is reduced toward a later stage equalizer. By way of 
example, FIG. 17 shows that the equalizer is divided in two 
parts where a prestage equalizer (multi-user equalizer) 71 
cancels an interfering component which is located outside an 
equalization range of a poststage, single user equalizer 21'. 
At this end, a preprocessing including a soft interference 
cancel and a linear filtering according to MMSE(minimum 
mean square error) criteria, for example, is made, and 
subsequently, the poststage equalizer 2T performs an equal- 
ization of a single user having a number of paths equal to Q. 

[0189] Even when an equalization takes place in a cascade 
manner and a linear filter is used in a prestage processing, it 
is possible to prevent the quantity of calculation from 
increasing prohibitively. 

[0190] An arrangement of a multiple output turbo -receiver 
according to an embodiment which is based on this funda- 
mental concept of the first aspect of the invention (2) of the 
turbo-reception method and an exemplary arrangement of 
MIMO system which incorporates the present invention is 
shown in FIG. 18. It is to be understood that parts corre- 
sponding to those shown in FIG. 1 are designated by like 
reference characters as used before, without repeating their 
description. (The same holds true in subsequent description.) 

[0191] Transmitted signals from each transmitter is 
received through transmission paths (channels) by a turbo- 
receiver 30. The received signal r(k) is input to a multi-user 
equalizer 71, from which signals from N transmitters are 
delivered in the form of signals u 1 (k), . . . , u N (k), each of 
which is provided in the form it is free from interferences by 
signals from other transmitters, and channel values a 1 (k), . 
. . , a N (k) to be input to single user equalizers 21-1, . . . , 
21-N. These SISO equalizers 21-1, . . . , 21-N deliver 
log-likelihood ratios A 1 [b 1 (k)], . . . , A-^b^k)]. Subsequent 
processing remains similar to FIG. 1. However, the channel 
values a. 1 (k), . . . , a. N (k) which are used in the single user 
equalizers 21-1, . . . , 21-N are channel values which are 
obtained subsequent to the multi-user equalization, and are 
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distinct from the channel matrix H. Accordingly, a 1 (k) 7 . . . 
, a N (k) are referred to as post-equalization channel infor- 
mation. 

[0192] The operation will now be described more specifi- 
cally. 

[0193] Equations (23) to (26) are defined in a similar 
manner as described above in connection with FIG. 1 in 
consideration of the number of multipaths (channels) Q. 

[0194] The purpose of the poststage equalizers 21-1, . . . 
,21-N shown in FIG. 18 is to equalize intersymbol interfer- 
ence channels by signal symbols of respective users them- 
selves [b n (k), b n (k-l), . . . , b n (K-Q+l)] (where n=l, . . . , 
N). At this end, the prestage equalizer 71 operates to 
eliminate interferences other than [b n (k), b n (k-l), . . . , 
b n (K-Q+l)] (where n=l, . . . , N) within y(k). This will be 
discussed below more quantitatively. 

[0195] Initially, using the a priori information X 2 p [b n (k)] 
(where n=l, . . . , N) of the equalizer 71 which is fed back 
from the decoders 24-1, . . . , 24-N, a soft decision trans- 
mitted symbol estimate b'(k) is determined according to the 
equation (15). 

[0196] This soft decision transmitted symbol b' n (k) and a 
channel matrix H are used to generate a replica HB'(k) of an 
interfering signal, which is then subtracted from the received 
matrix y(k). 



y'Jk)=y(k)-H-B'{k) 
=H-{B{k)-B'{k))+n{k) 

[0197] where 

BXk)=[b'^(k+Q-l) . . . b^(k) . . . b' T (k-Q+l)Y 

b'(k+q)=[b' 1 (k+q)b'^k+q) . . . b' n (k+q) . . . b'^(k+q)] 
^:q=Q-l, . . . , 1 

bXk+q)=[b' 1 (k+q)b' 2 (k+q) . . . 0 . . . b' ys (k+q)Y :q=0, . 
■ • , -Q+i 



(27) ' 

(28) ' 

(29) ' 
(53) 
(54) 



[0198] it being noted that b'(k+q) has a zero element at 
n-th position. 

[0199] The operation of subtracting the interference in this 
manner will be referred to hereafter as a soft interference 
cancel. Assuming that a replica of the interfering signal is 
generated in an ideal manner, it will be seen that y'n(k) 
resulting from the subtraction can have only the symbol 
b n (k) of the n-th user and an intersymbol interference 
component caused by the symbol [b n (k-l), . . . , b n (k-Q+l)] 
of the n-th user itself which results from putting the n-th 
element of b'(k+q) equal to 0 for q=l, . . . , -Q+l in the 
equation (54). 

[0200] In effect, a contribution from the signal of the n-th 
user (transmitter) to the received matrix r(k) is only that 
resulting that from the symbol [b n (k), b n (k-l), . . . , 
b n (k-Q+l)]. However, it will be understood from the defi- 
nition of the received matrix y(k) given by the equation (21) 
that a contribution from the signal of the n-th user (trans- 
mitter) within the received matrix y(k) which results as 
synthesis of the multipaths will contain, when based on the 
k-th symbol b n (k), intersymbol interference components 
caused by future symbols [b n (k+Q-l), b n (k+Q-2), . . . , 
b n (k+l)] . Thus, the interference replica includes interfer- 
ence components from the future symbols. In this respect, 



the difference matrix y'(k) defined by the equation (27)' is 
distinct from the difference matrix y'(k) defined by the 
equation (27). 

[0201] Accordingly, a next step in the prestage processing 
in the equalizer 71 is to eliminate the residue interference 
which remains after the soft interference cancel, namely, a 
residual interference component which results from an 
imperfect synthesis of the interference replica H-B'(k) and 
interference components between future symbols from 
y' n (k) using MMSE (minimum mean square error) criteria 
linear filter. In other words, this elimination takes place by 
an arrangement such that a filtering of y' n (k) using the filter 
characteristic w n (k) as indicated by the equation (55) is 
equal to a sum of the symbols [b n (k), b n (k-l), . . . , 
b n (K-Q+l)] each multiplied by channel values ct ln (k), 
a 2n (k), . . . , a Qn (k). 



>v n Ii (^-y' n (^-2^o Q - 1 a q (A:)-b n (^- ? )=a n H (A:)-6 n (A:) 



(55) 



[0202] Accordingly, what is required is to calculate the 
equation (55) by determining the filter characteristic w n (k) 
and the post-equalization channel value (channel informa- 
tion) cc n (k). The derivation of w n (k) and a n (k) will be 
described. It is to be noted that while the filter characteristic 
w n (k) is distinct from the filter coefficient w n (k) given by the 
equations (32) and (34), similar denotations will be used for 
purpose of convenience. 

[0203] Desired solutions are defined as solutions of the 
following optimization problem: 



(w n (k),a. n (k))= a ig min|K H (A)y n (fc)-a^*)|f 

[0204] provided a ln (k)=l. 



(56) 



[0205] In other words, w n (k) and a. n (k) which minimizes 
the right side of the equation (56) are determined. The 
constraint requirement a ln (k)=l is added in order to avoid 
solutions which result in cx n (k)=0 and w n (k)=0. While solu- 
tions can be obtained under the constraint requirement, 

K(*)|| 2 =i 

[0206] a solution for ct ln (k)=l will be described below 
herein. For brevity, the problem will be redefined. Namely, 
the right side of the equation (56) is defined as m n (k) which 
is minimized in terms of w and a. 



m n (£)=arg min||m n H (/:)-z n (&)|| 



(57) 



[0207] provided m n H (k)-e 

mq+i =— I (which is equivalent 

to ct ln (k)=l) and where 



^H^)-an(*)T 

z n (k^[ y J(k), K*) n T ] T 

e MQi-l = [0 . - . 1 . . . 0] T 

[0208] it being understood that e MQ+1 
(MQ+l-th position). 



(58) 
(59) 
(60) 

has "1" element at 



[0209] A solution of the optimization problem is given as 
follows according to Lagrange's method of indeterminate 
coefficients described in literature [2], S. Haykin, Adaptive 
Filter Theory, Prentice Hall, pp. 220-227; 



m n(^)-"^ZZ 1,e MQ+l/( e MQ+l H "^ZZ ^MQ+l) 



H. 



-1. 



[0210] where 



(61) 



(62) 
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[0211] E[A] representing an expected value (a mean 
value) of A. 



= E 



H ■ /\ n (k) ■ H M + a 2 1 H% 



(63) 



A n (£)=diag [Z) n (^+e-l) 5 . . . £>Jk) 9 . . . ^D n (A:-G + li$4) 

[0212] where I represents a unit matrix, and a 2 noise 
power (a covariance of white Gaussian noise). 



A![W]=log 



Pr[fr„(*) = +1 I u n (k),k=0, ... ,B] 
Pr[b n (k) = -1 | u n (k), k =0, ... ,5] 



(69) 



-X 1 [fe n (fc)]+^ p [i n (*)] 



(70) 



[0218] where X 2 [b n (k)] represents an extrinsic informa- 
tion fed to a succeeding decoder 24-n, and [b n (k)] a 
priori information applied to the equalizer 31. Decoder 24-n 
calculates the log-likelihood ratio A 2 according to the fol- 
lowing equation: 



fh^Q-i) 0 0 0 

hn(Q-2) hn(Q-l) 0 0 



h n (Q-l) 



(65) 



r[b n (i) = +1 | A x [*„(/)], i= 0, ... ,5] (71 ) 



Z) n (*+<7)=diag [l-6^(^) 3 
,l-// 2 N /r+<7]]:<7=2+l, ... ,1 

D n (k+q)=diag [l-b'\(k+q), . . 
:q=0, . . . 



,1, 



4-i' 2 N (^+^)] 



(66) 



(67) 



[0213] where diag represents a diagon al matrix (all ele- 
ments other than those located along the diagonal of the 
matrix being zeros). Thus, if the channel matrix H and a 2 are 
known, m n (k) can be determined according to the equation 
(61). Accordingly, w n (k) and ct n (k) can then be determined 
according to the equation (58). 

[0214] Using the filter characteristic w n (k) which is cal- 
culated in this manner, y' n (k) is filtered according to the 
following equation; 



(68) 



[0215] where represents a conjugate transposed matrix. 

[0216] These filtered n processed results are fed to corre- 
sponding equalizers 21-n which follow. In this manner, a 
received signal u n (k) which corresponds to the left side of 
the equation (1) from the n-th user is obtained, a mn (k) which 
corresponds to the channel value h mn (q) on the right side of 
the equation (1) is obtained, and the equation (55) which 
corresponds to the equation (1) is determined. Accordingly, 
cc n (k) is applied as an equalizer parameter (channel value) to 
a succeeding equalizer 21-n. This completes the prestage 
processing by the equalizer 71. 

[0217] A processing which takes place in the succeeding 
equalizer 21-n and thereafter will now be described. As 
mentioned previously, because the equation (55) corre- 
sponds to the equation (1), the operation which takes place 
in the equalizer 21-n for every user may proceed in the 
similar manner as that of the equalizer 21 shown in FIG. 31, 
and will not be repeated again since it is disclosed in 
literature [1], Each equalizer 21-n receives u n (k), a n (k) and 
a priori information K 2 [b n (k)] from the decoder 24-n as 
defined above, and calculates and delivers outputs a log- 
likelihood ratio A 1 (LLR) representing the ratio of the 
probability that each encoded bit is +1 and the probability 
that it is -1, as follows: 



Pr[b n (i)= -1 |Ai[Mi)l " = 0, ... , B] 

ss A 2 pb„(i)] + Af [MO] ( 72 ) 



=^n(0]+^/[^(0] 



(72) 



[0219] where X 2 [b n (i)] represents an extrinsic information 
applied to the equalizer 71 and the equalizer 21 during the 
iteration and the X ± F [b n (k)] a priori information applied to 
the decoder 24-n. The arrangement shown in FIG. 18 
performs an iterated equalization and decoding to improve 
an error rate. 

[0220] A functional arrangement of the multi-user equal- 
izer 71 will be briefly described with reference to FIG. 19. 
A received signal from each antenna is processed in a 
receiver 70 as a matrix r(k)=[r 2 (k) . . . r M (k)], which is used 
in a received matrix generator 311 to generate a received 
matrix y(k) according to the equation (21) which takes the 
individual multi-paths (channels) into consideration. 

[0221] On the other hand, the received signal r(k) from the 
receiver 70, and a known series signal such as a unique word 
series used for channel estimation which corresponds to 
each transmitter and which is fed from a unique word 
memory 29 are input to a channel estimator 28 in order to 
estimate a channel matrix H. 

[0222] A priori information >^ p [b n (i)], . . . ,^i P [b N (i)] are 
subtracted from the log-likelihood ratios A 2 [b-^i)], . . . 
,A 2 [b N (i)] delivered from the respective decoders 24-1, . . . 
, 24-N to derive extrinsic information ^ 2 [bi(k)], . . . , X 2 
[b N (k)] , which is then input to soft decision symbol esti- 
mators 313-1, . . . , 313-N where soft decision transmitted 
symbols b , 1 (k), . . . , b' N (k) are calculated according to the 
equation (15) and are then input to an interference matrix 
generator 72. In the interference matrix generator 72, a 
matrix B'(k) of symbol estimates which can be interference 
signals from other transmitters are generated for each n 
according to the equations (29)', (53) and (54). A product of 
these N matrixes B'(k) and the channel matrix H is generated 
by other interfering signal estimators 73-1, . . . , 73-N, 
respectively, thus determining the replica of interfering 
components H B(k). 

[0223] These N interfering component replicas H-B(k) are 
subtracted from the received matrix y(k) in subtractors 74-1, 
. . . , 74-N, respectively, thus providing difference matrixes 
y\(k), y' N (k). 

[0224] The soft decision transmitted symbols b , 1 (k), . . . 
,b' N (k) are input to an error matrix generator 75 where error 
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matrices A 2 (k), . . . , A N (k) are generated according to the 
equations (64), (66) and (67). These error matrices, the 
channel matrix H and the noise power o 2 are input to a filter 
characteristic estimator 76 where the filter characteristic 
w n (k) and the post-equalization channel information ot n are 
estimated according to the equations (58), (60), (61), (63) 
and (65). These filter characteristics w 1? . . . ,w N and 
difference matrixes y'^k), . . . , y' N (k) are multiplied together 
in filter processors 77-1, . . . , 77-N, respectively, or the 
difference matrixes are filtered, thus determining a compo- 
nent u 2 (k), . . . , u N (k) of the received signal for the symbol 
[b n (k), b n (k-l), . . . , b n (K-Q+l)] from each user and for 
each path, from which interferences from other user signals 
are eliminated. These components are fed, together with the 
post-equalization channel information cx 1 (k), . . . , ct N (k) 
which are determined in the filter characteristic estimator 76, 
to the single user equalizers 21-1, . . . , 21-N shown in FIG. 
18. 

[0225] A processing procedure for the turbo -reception 
method according to the first aspect of the invention (2) is 
shown in FIG. 20 where steps corresponding to those shown 
in the procedure of FIG. 3 are designated by like step 
numbers as used before. However, the calculation of the 
interference replica matrix B' n (k) which takes place at step 
S4 follows the equations (29)', (53) and (54). The step S13 
uses the soft decision transmitted symbol b' n (k) to generate 
the error matrix A n (k) according to the equations (64), (66) 
and (67). Step S14 uses the channel matrix H, the noise 
power a 2 and the error matrix A n (k) to determine the 
residual interference eliminating filter w n (k) and the channel 
information a n according to the equations (58), (60), (61), 
(63) and (65). Step S15 filters the difference matrix y' n (k) 
according to the residual interference eliminating filter char- 
acteristic w n (k) to determine u n (k). At step S16, a single user 
equalization is applied to each filtered result u n (k) to deter- 
mine the log-likelihood ratio A n [u n (k)], which is then 
decoded at step S10. In other respects, the procedure is 
similar to that shown in FIG. 3. 

[0226] In the forgoing description, the extent of equaliza- 
tion in the poststage equalizer 21-n is defined as a zone for 
intersymbol interference by the symbol [b n (k), b n (k-l), . . . 
, b n (K-Q+l)] (where n=l, . . . , N), but such extent of 
equalization is adjustable. For example, when Q has a very 
high value, a computational load on the poststage equalizer 
21-n will be much greater. In this instance, the extent of 
equalization by the poststage equalizer 21-n is chosen to be 
Q'<Q while the prestage equalizer 71 may be rearranged so 
as to eleminate intersymbol interferences between the signal 
of the same user outside the zone [b n (k), b n (k-l), . . . , 
b n (K-Q'+l)] (where Q'<Q and n=l, . . . , N). Such modifi- 
cation will be described later. When the equalization is 
divided into the prestage and poststage, a previous symbol 
memory 32 may be provided, as indicated in broken lines in 
FIG. 19, in association with the channel estimator 28 so that 
the hard decision transmitted symbol 6 n (k) may also be used 
in the estimation of the channel value, thus permitting the 
accuracy of estimation to be improved. 

[0227] In the example shown in FIG. 17, the prestage 
multiple output equalizer 71 equalizes transmitted signals in 
N series so that interferences from other series may be 
separated to provide signals u n of N series and the post- 
equalization channel information a n and subsequently, the 
signal u n of each of N series is processed by the post stage 



single user equalizer 22 -n so as to eliminate intersymbol 
interference of the same transmitted signal. In this manner, 
the equalization takes place in two stages which are in 
cascade connection. However, a cascade connection of three 
or more stages may also be used. 

[0228] By way of example, FIG. 21 shows that a received 
signal r m of M series with respect to a transmitted signal of 
N series is input to a first stage equalizer 81 to provide an 
equalized signal series er a (k) for 1st to U-th transmitted 
series from which interferencs by (U+l)th to N-th transmit- 
ted series is eliminated and its associated post-equalization 
channel information ecx(k), and an equalized signal series 
er 2 (k) for (U+l)-th to Nth transmitted series from which 
interferences by 1st to U-th transmitted series are eliminated 
and its associated post-equalization channel information 
ea 2 (k) while a second stage includes equalizers 82-1 and 
82-2. er 2 (k) and eot 1 (k) are input to the equalizer 82-1 where 
they are equalized to provide an equalized signal series 
er 3 (k) for 1st -U-L-th transmitted series among the 1st to Uth 
transmitted series from which interferences by (Uj+iyth to 
U-th transmitted signal are eliminated and its associated 
post-equalization channel information ea 3 , an equalized 
signal series er 4 (k) for (U 2 + l)-th to U 2 -th transmitted series 
among the 1st to U-th transmitted signal from which inter- 
ferences by the 1-st to U 1 -th transmitted series and U 2 -th to 
U-th transmitted signal are eliminated and its associated 
post-equalization channel information ea 4 (k), and an equal- 
ized signal series er 5 (k) for (U 2 +l)-th to U-th transmitted 
series among 1-st to U-th transmitted series from which 
interferences by 1-st to U 2 -th transmitted series are elimi- 
nated and its associated post-equalization channel informa- 
tion ect 5 (k). 

[0229] Similarly, the equalized signal series er 2 (k) and the 
channel information ea 2 (k) are input to the equalizer 82-2 in 
the second stage to provide an equalized signal series er 6 (k) 
and an associated post-equalization channel information 
ea 6 (k), and an equalized signal series er 7 (k) and its associ- 
ated post-equalization channel information eo. 7 (k). When 
N=5, equalizers 83-1 to 83-5 in a third stage represent single 
user equalizers shown in FIG. 18. Alternatively, an input 
equalized signal to the equalizer 83-3 may comprise two 
transmitted signals, and the equalizer 83-3 can eliminate 
mutual interferences between the two transmitted signals to 
provide a set of two equalized signals and their associated 
post-equalization channel information, which are in turn 
supplied to and equalized by next succeeding single user 
equalizers 84-1 and 84-2. As a further alternative, the 
equalizer 83-4 may receive the equalized signal er 6 (k) and 
the channel information ect 6 (k) to eliminate mutual interfer- 
ences for each of the transmitted signals which constituted 
the equalized signal er 6 (k), such as each of three transmitted 
signals, each of which may be interfered by two other 
transmitted signals as well as an intersymbol interference of 
each transmitted signal due to multipaths. As another alter- 
native, one or more of the equalizers 82-1 and 82-2 in the 
second stage may be arranged so that an equalized signal 
may be simultaneously obtained for each of a plurality of 
transmitted signals. 

[0230] What has been described above can be generalized 
by stating that equalizers in a first stage deliver a plurality of 
equalized signal series and a set of post-equalization channel 
information, and one or more equalizers may be provided in 
each of one or a plurality of stages which are in cascading 
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connection for each equalized signal series and its associated 
set of post-equalization channel information so that an 
equalized output or a log- likelihood ratio A 2 [b n (k)] may be 
delivered finally for each of the 1st to N-th transmitted 
series. 

[0231] When the equalization takes place in multiple 
stages which are in cascade connection, it is preferred that 
the number of paths Q for which interferences are to be 
eliminated be reduced toward a later stage so that the 
quantity of calculation can be reduced. In this instance, an 
interferring component from a path which disappears in a 
later stage be eliminated in an immediately preceding equal- 
izer stage. 

[0232] An equalization processing which occurs in the 
arrangement of FIG. 21 when a first stage equalizer 81 deals 
with N transmitted signals, each of which has a number of 
multi-paths equal to Q to provide a group of equalized signal 
series er 1 (k) comprising U transmitted signals and their 
associated post-equalization channel information ea. 1 (k) and 
a later stage equalizer 82-1 performs an equalization for a 
number of multipaths equal to Q' for each transmitted series 
will be described. 

[0233] In a similar manner as described above in connec- 
tion with the embodiment shown in FIGS. 18 and 19, an 
interference matrix generator 72 generates an interference 
matrix B'(k), but the equations used are changed from the 
equations (53) and (54) to equations (53), (54)' and (73) as 
follows: 



bXk+q)=[b' ± (k+q)b' 2 (k+q) . . . b'Jk+q) . . . b'^k+qft 
^:q=Q-\, ... 1 

b%k+q)=[0 ... 0 b'v +1 {k+q) . . . b'^(k+q)Y :q=0, . . . 

-G'+i 



bXk+q)=[b\(k+q)b'^k+q) 



b'Jk+q) . . . b' N (k+qy] 



(53) 
(54)' 



(73) 



[0234] The equation (54)' is intended to provide symbols 
for a first to U-th transmitted series themselves, and to 
provide an equalization except for intersymbol interference 
caused upon each series by itself and relative to each other 
due to multipaths Q' while the equation (73) is intended to 
eliminate the intersymbol interferences upon the 1st to U-th 
transmitted series themselves and relative to each other due 
to (Q'+l)-th to Q-th path in as much as the number of 
multipaths is reduced to Q f in the poststage equalization. 

[0235] The interferences matrix B'(k) which is obtained in 
this manner is used to generate an interference signal replica 
H B'(k) which is then subtracted from a received matrix y(k) 
as follows: 

y ' S (k)=y(k)-H-B'(k) (27)" 
=H- {B{k)-B '(£))+«(£) (28) " 

[0236] This operation of subtracting the interference is 
referred to hereafter as a soft interference cancel. Assuming 
that a replica HB'(k) of an interfering signal is generated in 
an ideal manner, it will be seen that y g (k) can only have 
signal components for symbols [b n (k), b n (k-l), . . . , 
b n (k-Q'+l)] (where n=l to U) for 1st to U-th transmitted 
symbols. 

[0237] A residue of interferences which remain after the 
soft interference cancel is eliminated with a linear filter of a 
MMSE criteria in a similar manner as mentioned previously. 
In this instance, the equation (55) is replaced by an equation 
(55)' indicated below. 



^(^■> ? ' g ( A: )- 2 n=i U 2 q .o Q '- 1 « nq (^)-b 11 (^-^)= 



(55)' 



[0238] where 



a s (k)=[o. uo (k), . . . rO^v-ifc)* ■ - ■ > a u,o(#) 5 ■ ■ ■ 

) a U,Q'-lW] 

(?+i)T 



(55-1) 



(55-2) 



[0239] The derivation of w g (k) and ct g (k) takes place in a 
similar manner as described previously to determine w g (k) 
and oc g (k) which minimizes the right side of the following 
equation, which stands for the equation (56): 



(>v g (A:),a g (*))=arg min\\w™(k)y> g (k)-ai™ (*)-fc 8 (*)f 

[0240] provided a ± 0 (k)=l. 



(56)' 



[0241] The constraint requirement is added in order to 
avoid solutions which may result in a g (k)=0 and w g (k)=0. 
While a constraint requirement that||a g (k) || 2 =1 may also be 
used to provide a solution, in the description to follow, the 
problem is rewritten as follows for o. 10 (k)=l: 



[0242] ti m g (*)=arg min||m g H (/:)-z g (^)|| : 
[0243] provided m g H (k)-e MQ , +1 =-l 
[0244] where 



(57)' 



m 8 (A)4 Wg T (*)-a g T (*)] 

* g (*>b/(*)A*)gT 



-MQ'+rt^ . . . 1 . . . Of 7 

[0245] it being understood that e MQ , +1 
(MQ'+l)-th position. 



(58) ' 

(59) ' 

(60) ' 

has "1" element at 



[0246] A solution of the optimization problem is given as 
follows according to Lagrange's method of indeterminate 
coefficients disclosed in literature [2]; 



m g (k) — R zz '^MQ'+l/^MQ'+l '^ZZ " e MQ'+l) 

[0247] where 

R zz =E[z g (kyz*(w 

[0248] 

R zz = E[Z g (k)-^(k)] 



(61)' 



(62)' 



= E\ 



(62 y 

(63 )' 



U*)=diag[2) n (/:+<2-l), . . . ,D n (*), . . . J> n (*H2+l)l64)' 
[0249] 



«Ke-u 

hnQ-3) 



0 



i h 



1(0) 



0 
0 

^1(0-1) 

■ 



hlJ(Q-2) 
hu(Q-3) 

■ 



0 



0 
0 

m 

hu(Q'-l) J 



(65 y 



D^k+q^dmgll-b^^k+q), . . . f l-b a n (k+q) 9 . . . 
,l-b%(k+q)];q=Q + l, ... ,1 

D n (*+?)=diag[l, . . . ,l,l-b% +1 (k+q), . . . ,l-b%(k+ 
q)l-q=0, . . . -Q'+l 

D n (k+q)=di a g[l-b' 2 1 (k+q), . . . ,l-b'\(k+q), . . . 
,l-b%(k+q)l-q=Q' f . . . -Q+l 



(66) 

(67)' 
(74) 
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[0250] Thus when channel parameters are known, m g (k) 
can be determined according to the equation (61)'. In addi- 
tion, w g (k) and ct g (k) (=ect 1 (k)) can be determined accord- 
ing to the equation (58)\ Such calculations can be made in 
the filter characteristic estimator (76) shown in FIG. 19, for 
example, and the following equation is calculated in the 
filter processor 77-1 for purpose of filtering: 

[0251] This equalized output er 2 (k) and the post-equaliza- 
tion channel information ea. 1 (k)=a g (k) are fed to the post- 
stage equalizer 82-1. 

[0252] When there are five transmitted series (users), for 
example, which are divided into a group of three transmitted 
series (users) and a group of two transmitted series (users), 
in the manner mentioned above, the above algorithm is 
carried out with U=3 and 2, and the two equalized outputs 
er 2 (k), ect 1 (k) and er 2 (k), ect 2 (k) are input to equalizers 
which are designed to deal with the three transmitted series 
and the two transmitted series, respectively, thus obtaining 
an equalized output for each transmitted series. 

[0253] Reflecting an error correction decoding result for a 
signal being detected into a soft decision transmitted symbol 
in the manner mentioned above is also applicable to a single 
user turbo-equalizer receiver shown in FIG. 8, RAKE 
synthesis turbo -receiver shown in FIG. 9, a turbo -receiver 
including an adaptive array antenna receiver shown in FIG. 
10, a generalized turbo -receiver including a channel esti- 
mator 42 shown in FIG. 12. 

[0254] In FIGS. 13, 14 and 15, a symbol hard decision 
value which has been determined as being as likely to be 
certain is also used as a reference signal in the estimation of 
the channel matrix H and the covariance matrix U during a 
second and a subsequent iteration. However, during a second 
and a subsequent iteration, only the unique word may be 
used as a reference signal to utilize the equation (51) to 
estimate the covariance matrix U while omitting the esti- 
mation of the channel estimation and the estimation of the 
covariance matrix U which utilize the symbol hard decision 
value. 

[0255] First Aspect of the Invention (2) (Parallel Trans- 
mission) 

[0256] There is a proposal that information series c(i) from 
a single user be transmitted in a plurality of parallel series in 
order to achieve a high rate transmission with a high 
frequency utilization efficiency. An embodiment of a turbo- 
receiver incorporating the present invention which may be 
used for such transmitted signal will now be described. 

[0257] Referring to FIG. 22, where parts corresponding to 
those shown in FIG. 1 are designated by like reference 
characters as used before, on the transmitting side, a modu- 
lated output signal b(j) from a modulator 13 is input to a 
series-parallel converter 14 where each symbol b(j) is 
sequentially distributed into N series. It is assumed that there 
are series signal b 1 (k), . . . , b N (k), the number N of which 
is equal to or greater than 2. While not shown, these signals 
are transmitted from N antennas after convertion into radio 
frequency signals. 

[0258] These N series signals are propagated through 
channels (transmission paths) to be received by the turbo- 
receiver according to the present invention. The receiver has 



one or more receiving antenna, and the received signal is 
input to a multiple output equalizer 31 as a baseband digital 
received signal r m (k) (where m=l, 2, . . . , M) including one 
or more (M) signals. The received signal r m (k) is generated 
in a manner as shown in FIG. 30B, for example. 

[0259] The multiple output equalizer 31 is constructed in 
the same manner as shown in FIG. 2, and performs accord- 
ing to a processing procedure as shown in FIG. 3. Accord- 
ingly, an extrinsic information ^i[b(i)] is subtracted from a 
log -likelihood ratio A 2 [b(i)] from a decoder 24 shown in 
FIG. 22 in a subtractor 25, and the subtracted output is 
interleaved by an interleaver 26 to provide a priori infor- 
mation ^ 2 [b(j)], which is then converted in a series-parallel 
converter 15 into N series of a priori information X, 2 [b 2 (k)], 
. . . , X 2 [b N (k)] to be input to the multiple output equalizer 
31. 

[0260] Accordingly, N series of received signals are sub- 
ject to a linear equalization in the multiple output equalizer 
31 in the similar manner as mentioned previously, delivering 
N log-likelihood ratio series A ± [b a (k)], . . . , A 1 [b N (k)], 
which are then input to a parallel-series converter 16 to be 
converted into a single log-likelihood ratio series A-^b^)] to 
be supplied to a subtractor 22. With this arrangement, the 
input signal format to the multiple output equalizer 31 is 
similar to that described in connection with FIGS. 1 to 3, and 
accordingly, N series log-likelihood ratios A 1 [b 1 (k)], . . . , 
Aj_ [b N (k)] can be obtained by the equalization which has 
been mentioned above with reference to FIGS. 1 to 3, and it 
will be readily seen that an iterative decoding processing is 
permitted by the use of the series-parallel converter 15 and 
the parallel- series converter 16. In a manner corresponding 
to a transmitted signal from a n-th transmitter in the arrange- 
ment of FIGS. 1 to 3, an n-th or (n-th column) transmitted 
signal among N parallel transmitted signals will be equal- 
ized. It will be readily seen that the embodiment described 
above in connection with FIGS. 4 to 7 is applicable to the 
reception of the parallel transmission of the N series signals. 
By processing in a plurality of equalizer stages which are in 
cascade connection as illustrated in FIGS. 18 to 21, the 
reception characteristic can be improved as compared to the 
processing through a single equalizer stage as shown in 
FIGS. 1 to 31. 

[0261] The turbo -reception method and the turbo-receiver 
according to the present invention are also applicable to the 
reception of convoluted code/turbo-code+interleaver+multi- 
value modulation such as QPSK, 8PSK, 16QAM, 64QAM 
etc., TCM (trellis coded modulation)/turbo TCM. 

[0262] Generation of M Received Signals 

[0263] M received signals r 2 (k), . . . , r M (k) are derived 
from M antennas #1, . . . , #M, but may be derived from a 
single antenna. Alternatively, M (which is greater than L) 
received signals may be obtained from L (which is an integer 
equal to or greater than 2) antennas. While not specifically 
shown in FIG. 1, received signals from antennas #1, . . . , 
#M are converted into baseband received signals r 1 . . . , r M 
in a baseband converter and sampled to provide digital 
signals r 1 (k), . . . , r M (k) at discrete time k. 

[0264] As shown in FIG. 30£>, for example, received 
signals received by L=2 antennas #1 and #2 may be con- 
verted into baseband signals in baseband converters 61-1 
and 61-2, respectively, the outputs of which may be sampled 
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by A/D converters 64-1, 64-2 and 64-3, 64-4 in accordance 
with a sampling signal from a sampling signal generator 62 
and a sampling signal which is obtained by shifting the 
phase of the first mentioned sampling signal by one half 
period thereof T/2 in a phase shifter 63 to provide digital 
signals r 1 (k), r 2 (k), r 3 (k), r 4 (k), respectively. Then the digital 
signal may be input to a turbo-receiver 30 as shown in FIG. 
1 or FIG. 18 or FIG. 22 to provide N decoded outputs. It 
should be understood that the frequency of the sampling 
signal from the sampling signal generator 62 is chosen so 
that each sampling period of the received signals r 1 (k), . . . 
, r 4 (k) input to the turbo -receiver 30 coincides with the 
sampling period which is used when a single received signal 
r M (k) is received per antenna. 

[0265] Effects of the Invention 

[0266] As discussed above, according to the first aspect of 
the invention (1), there is realized a multiple output (MIMO) 
reception method. To illustrate a quantitative effect, an error 
rate response is graphically shown in FIGS. 23 and 24. In 
each Figure, E b /N 0 on the abscissa represents a bit power to 
noise ratio. A simulation presumes following conditions: 



Number of users (transmitters) 

The number of multipaths of each user Q 

Number of reception antennas 

Number of information symbols per frame 

Number of unique words per frame 

Channel estimation method 

Error correction code 

Doppler frequency 
Modulation 
Transmission rate 
Decoder 24 
Number of iterations 
No fading within a frame 



2 
5 
2 

450 bits 
25 bits 

RLS (obliteration factor 
0.99) 

rate Vz } constrained length 3 

convoluted code 

1000 Hz (Rayleigh fading) 

BPSK 

20 Mbps 

Max- Log- Map decoder 
4 



[0267] No approximation by the matrix inversion lemma 
for the inverse matrix is used in the calculation of the filter 
coefficient w. 

[0268] FIG. 23 shows an error rate characteristic when the 
channel estimation is achieve perfectly (without no estima- 
tion error) or when the channels are known for a number of 
users N=2, a number of reception antennas M=2, a number 
of Raleigh paths Q=5. It will be seen that initial run is not 
an iteration and the second iteration is the first iteration. It 
will be seen that the error rate characteristic is significantly 
improved by the iteration. It will be seen from this that the 
turbo -reception method according to the invention for 
MIMO operates in a proper manner. 

[0269] FIG. 24 illustrates the effect of an iterative channel 
estimation (according to the fourth aspect of the invention). 
The abscissa indicates a threshold value Th. E b /N 0 is fixed 
to 4 dB (where E b is for one user), and Th=1.0 may be 
considered as a conventional method in which no symbol 
hard decision value is selected, or where a channel estima- 
tion using a symbol hard decision value is not made. In this 
instance, it will be seen from the drawing that BER char- 
acteristic has a little effect of iteration because the channel 
estimation is inaccurate. A threshold Th=0 indicates that all 
of hard decision values are directly used, and when hard 



decision values of information symbols are used in this 
manner, it will be apparent from the drawing that a mean bit 
error rate is improved, allowing the channel estimation to be 
performed exactly to a corresponding degree. For a thresh- 
old value on the order of 0.2 to 0.6, the mean bit error rate 
is reduced than when Th=0,indicating that it is preferable 
that only a hard decision value which is likely to be certain 
be used. It will be seen that Th which is around 0.25 is most 
preferable. 

[0270] FIG. 25 shows an error rate characteristic of an 
MIMO reception method in which a transmitted symbol 
hard decision value which is determined to be likely to be 
certain according to a threshold value is used in the channel 
estimation or which employs an iterative channel estimation 
in a form of a curve 66. In this instance, the threshold is set 
up at 0.25 and the result shown represents a characteristic 
after four iterations with parameters that N=2, M=2, 
Q=5Rayleigh, f d T s = 1/20,000, and 900 symbols/frame. For 
purpose of comparison, an error rate characteristic with a 
perfect channel estimation is shown as a curve 67 while an 
error rate characteristic in which a hard decision value of 
information symbol is not used in a channel estimation or 
only a single channel estimation is made without iteration is 
shown by a curve 68. It will be seen from this graphic 
illustration that when the iterative channel estimation is 
used, the error rate characteristic approaches to that obtained 
by a perfect channel estimation. 

[0271] With the channel estimation method mentioned 
above, by determining whether a hard decision value is or is 
not likely to be certain on the basis of a decoded soft 
decision value, and by using symbol information having a 
hard decision value which is likely to be certain in the 
channel estimation during the next iteration, the channel 
estimation can be performed more correctly, allowing a 
decoding quality to be improved. 

[0272] In order to confirm the effect of an embodiment in 
which a covariance matrix U (for noise other than Gaussian 
noise) is estimated, a simulation is made with parameters 
indicated below. 



Number of users (transmitters) 

The number of multipaths of each user Q 
Number of reception antennas 
Number of information symbols per frame 
Error correction code 

Doppler frequency 
Modulation 
Transmission rate 
Decoder 24 
Number of iterations 



3 (one of which is chosen as 
an unknown interference: i(k)) 
5 
3 

450 bits 

rate Vi, constrained length 3 

convoluted code 

1000 Hz (Rayleigh fading) 

BPSK 

20 Mbps 

Log-MAP decoder 
4 



[0273] Three users (transmitters) are chosen to be of an 
equal power. 

[0274] FIG. 26 shows a result of simulation of BER (bit 
error rate) characteristic of a turbo-receiver which estimates 
H and U as shown in FIGS. 14, 15 and 16, and FIG. 27 
shows BER characteristic which directly uses the turbo- 
receiver shown in FIG. 1 (which employs the method shown 
in FIG. 13). In FIG. 26, it is assumed that the noise 
comprises only white Gaussian noise, and it is seen that two 
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or more iterations of the channel estimation and the decod- 
ing processing has little effect. By contrast, it will be seen 
from FIG. 27 that as the number of iterations is increased, 
an improvement in the BER characteristic is achieved, and 
in addition, BER is considerably less than that shown in 
FIG. 26 for the same E b /N 0 . 

[0275] In order to confirm the effect of the embodiment 
(according to the second aspect of the invention) in which an 
error correction decoding result is reflected into a symbol 
soft decision value b' n (k) of a received signal from an 
intended user (transmitter), a simulation is made with 
parameters as indicated below 



Number of users (transmitters) N 4 

The number of multipaths of each user Q 5 

Number of reception antennas 2 

Number of information symbols per frame 900 bits 

Error correction code convoluted code (coded rate: 

L /2, constrained length 3) 

Modulation BPSK 

Decoder 24 Log-Map decoder 

Erroneous coded rate V2 

Number of iterations 5 
A choice is made that f(b' n (k))= axb' n (k) 



of receiving antennas or the like. In addition, a different 
value may be chosen as an optimum value for ol. 

[0277] When BPSK modulation is used for a number of 
users (transmitters) which is equal to N, a number of 
multipaths from each transmitter which is equal to Q and a 
number of receiver antennas which is equal to M, the 
quantity of calculation which is required in an equalizer 
when a conventional single user turbo -receiver is directly 
extended to a multiple output (MIMO) is on the order of 
2 N(Q_1) , as mentioned previously, but with the turbo-recep- 
tion method according to the third aspect of the invention, 
the quantity of calculation can be reduced to the oreder of 
N(MQ) 3 . By way of example, assuming that N=8, Q=20 and 
M=8, 2 N(Q - 1) ~5-10 45 while N(MQ) 3 ~3710 7 , thus demon- 
strating that the quantity of calculation can drastically be 
reduced according to the turbo-reception method according 
to the second aspect of the present invention. 

[0278] A simulation has been conducted under conditions 
given below in order to confirm that a good bit error rate 
characteristic can be obtained according to the turbo-recep- 
tion method according to the third aspect of the invention. It 
is assumed that the channel matrix H is known. 



[0276] FIG. 28 shows BER characteristic of a multiple 
output turbo-receiver shown in FIG. 1 and a multiple input 
multiple output turbo -receiver in which an error correction 
decoding result is reflected into b' n (k), with plotted points 
shown in black for the former and in white for the latter. It 
is to be noted that in the plot, a circle represents an initial 
run, a down-directed triangle a second iteration, a diamond 
a third iteration, a left-directed triangle a fourth iteration and 
a right -directed triangle a fifth iteration. FIG. 28 A shows a 
result of a simulation for the BER characteristic plotted 
against E b /N Q when a is fixed to 0.2, and FIG. 28B shows 
a result of simulation for the BER characteristic plotted 
against a when E b /N 0 is 6dB. It will be seen that a=0 
represents that b' n (k)=0. It is noted from FIG. 28A that in the 
multiple input multiple output receiver in which an error 
correction decoding result is reflected into b' n (k), an 
improvement in BER after a third and a subsequent iteration 
is greater in comparison to BER obtained during a preceding 
iteration, as compared with the multiple input multiple 
output turbo-receiver shown in FIG. 1, and that after a third 
and a subsequent iteration, when E^/Nq required to achieve 
each BER in a range BER>10~ 4 is compared, the multiple 
input multiple output turbo -receiver in which an error cor- 
rection decoding result is reflected into b' n (k) has a gain of 
about 0.5 dB or greater in comparison to the multiple input 
multiple output turbo -receiver shown in FIG. 1. It is also 
noted that during the fifth iteration at E b /N Q =6dB, BER= 
10~ 5 is attained, which is a reduction in BER by a factor of 
1/10 or more in comparison to the receiver shown in FIG. 
1. It will be seen from FIG. 28B that an improvement is 
achieved in a range of a. indicated by an equality 0<a<0.6 
and that when a exceeds 0.6, the BER characteristic 
becomes degraded, preventing a correct decoding result 
from being obtained. From this result, it is seen that an 
optimum value of a is 0.2. However, it should be understood 
that the value of a. is not limited to the optimum value, but 
a suitable range of a. having an improving effect can be 
changed depending on the number of users to be received, 
a propagation environment having interferences, a number 



Number of users N 

The number of multipaths of each user Q 
Number of reception antennas 
Number of information symbols per frame 
Error correction code 

Doppler Frequency 
Modulation 
Decoder 24 
Transmission rate 
Decoder 

Number of iterations 



4 

5 
2 

900 bits 

rate: l /2, constrained length 3 

convoluted code 

1000 Hz (Rayleigh fading) 

BPSK 

Log-Map decoder 
20 Mbps 

Log-MAP decoder 
6 



[0279] It is assumed that the channel estimation takes 
place in an ideal manner. 

[0280] FIG. 29 shows a result of simulation for BER (bit 
error rate characteristic). The abscissa represents mean E b 
(bit power)/ N 0 (noise power), denotation fd in the graph 
represents a doppler frequency and D s the period of a trans- 
mitted symbol. MRC shown on this graph is BER charac- 
teristic upon Viterbi decoding of a signal after a maximal 
ratio combining on an order 10 (two antennasx5 paths) 
diversity channels, and this corresponds to the BER char- 
acteristic which results when the equalizer has perfectly 
canceled interferences. Thus, the quality of the receiver can 
be evaluated by seeing how the BER after the iterations is 
located close to the MRC curve. It is seen from FIG. 29 that 
with the turbo -reception method according to the second 
aspect of the present invention, the higher E b /N 0 , the less 
BER, and the greater the number of iterations, the more the 
BER characteristic approaches MRC BRE characteristic, in 
particular BER characteristic is located very close to MRC 
with a number of iterations equal to 6. Thus, it is confirmed 
that the multiple output turbo- receiver which employs the 
turbo-reception method according to the third aspect of the 
present invention operates properly under severe conditions 
of four users, each having five paths and using two reception 
antennas. 
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What is claimed is 

1. A turbo-reception method of receiving a signal from N 
transmitters where N is an integer equal to or greater than 2, 
comprising the steps of 

calculating a channel value h mn (q) and a channel matrix 
H from M received signals r m (where M is an integer 
equal to or greater than 1) and a known signal, where 
m=l, . . . , M; n=l, . . . , N and q=0, . . . , Q-l, and Q 
represents a number of multipaths of each transmitted 
wave; 

determining a soft decision transmitted symbol b' n (k) 
from N a priori information X 2 [b n (k)] where k repre- 
sents a discrete time; 

using the channel value h mn (q) and the soft decision 
transmitted symbol b' n (k), calculating an interference 
component HB'(k) upon a transmitted signal from an 
n-th transmitter, 

where 



H = 



\ H (0) ■■■ H(Q-l) 0 
L 0 W(0) ■■■ H[Q-D 

r /ill (q) ■■■ hlN{q) 1 



H (q) = 



' h 



MI(q) 



h>MN(q) J 



B'ikyW^ik+Q-l) . . . b'^k) . . . b' T (k-Q+l)Y 

b'{k + q)=[b\{k+q)b'^k+q) . . . b^(k+q)f 
q=Q-l . . . -Q+l q*Q 

b'(k)=[b\(k) . . .o . . . b' N (k)Y q=o 

b'(k) having a zero element at an n-th position and [] T 
representing a transposed matrix; 

subtracting the interferenced component H-B'(k) from a 
received matrix y(k) to obtain a difference matrix y'(k) 
where 



y(k)=[rT(k + Q-l)rT(k + Q-2) . 
r(k)=[r^k)r^k) . . . r M (k)Y 



using the channel matrix H or a reference signal to 
determine an adaptive filter coefficient w n (k) which is 
to be applied to a received signal corresponding to a 
transmitted signal from an n-th transmitter in order to 
eliminate residual interference components in the dif- 
ference matrix y'(k); 

and filtering the difference matrix y'(k) with the adaptive 
filter coefficient w n (k) to provide a log- likelihood ratio 
for the received signal which corresponds to the trans- 
mitted signal from the n-th transmitter and from which 
interferences are eliminated. 
2. A turbo -reception method according to claim 1 in which 
denoting a covariance matrix of noise components in the 
received matrix y (k) by U, using the soft decision trans- 
mitted symbol b' n (k) and the channel matrix H to calculate 
the adaptive filter coefficient w n (k) according to the follow- 
ing formula: 

wJk)=(HG (£)f/ H + U)' x h 

G(k)=diag[D(k+Q-l) . . . D(k) . . . D(k-Q+1)] 



Dik+q^dmgll-b'i^k+q), 
,l-b'^(k + q)] 

q=Q-l . . . -Q+l,q*0 

=dmg[l-b' 2 1 (k+q), ... ,1, 



, l-b'\(k+q), 



,l-b%(k+q)]q=0 



h = 



Hl,(Q-l)N+n "I 
i 

Hm ,{Q-l)-N+n J 



where H 1( - Q _ 1) . Nr+I1 represents an element of the matrix H 
located at row 1 and column (Q-l)N+n. 

3. A turbo-reception method of receiving a signal from N 
transmitters where N is an integer equal to or greater than 2, 
comprising the steps of 

calculating a channel value h mn (q) and a channel matrix 
H from M received signals r m where M is an integer 
equal to or greater than 1 and where m is 1, . . . , M; 
n=l, . . . . N; q=0, . . . , Q-l, and Q represents a number 
of multipaths of each transmitted wave; 

determining a soft decision transmitted symbol b' n (k) 
from N a priori information X 2 [b n (k)] where k repre- 
sents a discrete time; 

using the channel value h mn (q) and the soft decision 
transmitted symbol b' n (k) to calculate an interference 
component H-B'(k) upon a transmitted signal from the 
n-th transmitter where 



// = 



H{q) - ! 



H {Q) ... H (Q - U 
L 0 ^(0) 

•■ h\.N{q) 
■ h>MN(q) 



0 



H (Q-D J 



f ^11 {q) 



B\k)=[b fY {k+Q-l) . . . i>' T (Ar) . . . £' T (£-<2+l)J r 
bXk+q)=[b\(k+q)b' 2 (k+q) . . . b'^(k+q)Y 
q=Q-l . . . -Q+l q*Q 

bXk)=lb\(k) . . . -f(b' n (k)) ...b' N (k)Y q=0 

b'(k) having an element f(b' n (k)) at an n-th position, f( ) 
being a function of a variable b' n (k) and satisfying that 
f(0)=0 and 

d{f(b' n (k))}/d{b' n (k)}^0, and [] T representing a trans- 
posed matrix; 

subtracting the interference component H-B'(k) from a 
received matrix y(k) to determine a difference matrix 
y'(k) where 

yi^lr^k+Q-iy^k+Q-T) . . . r T (k)Y 

K*H'i(*> 2 (*) ■ ■ ■ >* M (*)] T 

using the channel matrix H or reference signal to deter- 
mine an adaptive filter coefficient w n (k) which is to be 
applied to a received signal corresponding to the trans- 
mitted signal from the n-th transmitter in order to 
eliminate residual interference components which 
remain in the difference matrix y'(k); 
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and filtering the difference matrix y'(k) with the adaptive 
filter coefficient w n (k) to provide a log- likelihood ratio 
for the received signal which corresponds to the trans- 
mitted signal from the n-th transmitter and from which 
interferences are eliminated. 

4. A turbo-reception method according to claim 3, further 
comprising, denoting a co variance matrix of noise compo- 
nents in the received matrix y(k) by U, using the soft 
decision transmitted symbol b' n (k) and a channel matrix H to 
calculate the adaptive filter coefficient w n (k) according to 
the following formula: 

G(&)=dkg[Z>(jfc+0-l) . . . D{k) . . . D(k-Q+1)] 

ZX^)=diag[l-&'\(&+<7), - - - , l-b'\{k+q), . . . 
,l-b%(k+£>] 

q=Q-l . . . -Q+l,q*0 

=diag[l-b a ± (k+q), . . . > l-b' 2 tl _ 1 (k%l+2E\f(b' n (k)']+E 
l-b'\ +1 {k), . . . ,l-b%(k+q)]q=0 

where E[] represents a mean value, and 



h = 



! H2,(Q-l)-N+n ! 

I : ! 

i ■ i 

i i 

L HMAQ-iyN+n J 



where H 1 f ( Q _i). N+n represents an element of the matrix H 
at row 1 and a column (Q-l) N+n. 

5. A turbo-reception method according to claim 2 or 4 in 
which an inverse matrix calculation during the calculation of 
the adaptive filter coefficient w n (k) takes place by using 
matrix inversion lemma. 

6. A turbo-reception method according to claim 1 or 2 in 
which the co variance matrix U of noise components in the 
received matrix y (k) is defined as o 2 I which is determined 
by a variance o 2 of a Gaussian distribution and a unit matrix. 

7. A turbo-reception method according to claim 1 or 2 in 
which the co variance matrix U of noise components in the 
received matrix y (k) is derived, using the received matrix 
y(k) and the estimated channel matrix H, as follows: 

f>=2 k .o Tr 0(^)-^(^)-0(^)-^^(A)) H 
B(k)=[b T (k+Q-l) . . . b T (k) . . . b T (k-Q+l)Y 
b(k+q)=[b 1 (k+q) . . . b lss (k+q)Y 
(q=-Q+l . . . Q-l) 

where Tr represents the length of the reference signal. 

8. A turbo -reception method according claim 2 or 3 in 
which D(k+q) is approximated by 0 for q^0, and D(k+q) is 
approximated by diag [0, . . . ,1, . . . , 0] for q=0. 

9. A turbo-reception method of receiving a signal from N 
transmitters where N is an integer equal to or greater than 2, 
comprising the steps of 

calculating a channel value h mn (q) and a channel matrix 
H from M received signals r m (where M is an integer 
equal to or greater than 1) and a known signal where 
m=l, . . . , M; n=l, . . . , N; q=0, . . . , Q-l, and Q 
represents a number of multipaths of each transmitted 
wave; 



determining a soft decision transmitted symbol b' n (k) 
from N a priori information X 2 [b n (k)] where k repre- 
sents a discrete time; 

using the channel value h mn (q) and the soft decision 
transmitted symbol b' n (k) to calculate an interference 
component HB'(k) upon a transmitted signal from an 
n-th transmitter where 



(0) 



... Htn- 



(2-1) 



0 



H = ! 



L 0 ^(0) ■■■ H{Q-L) J 

i Ml (q) ■■- hlN{q) 1 

H {q) = ! ; ■-. ; 

!_ tlMHq) ■■■ hMN{q) J 

B'(k)=\b' T (k+Q-X) . . . b' T (k) . . . b' T (k-Q+l)Y 
b'(k+q)=[b , 1 (k+q)b f 2 (k+q) . . . b'^(k+q)Y 
q=Q-l . . . -Q+l q*Q 

b'{k)=\b\{k) . . . o . . . b'^(k)Y q=o 

b'(k) having a zero element at an n-th position and 
[^representing a transposed matrix; 

subtracting the interference component H-B'(k) from a 
received matrix y(k) to determine a difference matrix 
y'(k) where 

y(k)=[^(k + Q-l)r^(k + Q-2) . . . r^(k)Y 
K*H'i(*> 2 (*) ■ ■ ■ r M (k)Y 

defining a covariance matrix of noise components in the 
received matrix y(k) as a 2 I which is determined from 
a variance a 2 of a Gaussian distribution and a unit 
matrix I, and filtering the difference matrix y'(k) with an 
adaptive filter coefficient w n which is determined by 



h = 



filAQ-l)-N-m 

H 2 ,(Q-l).N+n 



H M -QAQ-].).N+n 



to provide a log-likelihood ratio for the received signal 
which corresponds to the transmitted signal from the 
n-th transmitter and from which interferences are elimi- 
nated. 

10. A turbo -reception method of receiving a signal from 
N transmitters where N is an integer equal to or greater than 
2, comprising the steps of 

calculating a channel value h mn (q) and a channel matrix 
H from M received signals r m (where M is an integer 
equal to or greater than 1) and a known signal where 
m=l, . . . , M; n=l, . . . , N; q=0, . . . , Q-l, and Q 
represents a number of multi paths of each transmitted 
wave; 

determining a soft decision transmitted symbol b' n (k) 
from N a priori information X 2 [b n (k)] where k repre- 
sents a discrete time; 

calculating an interference component H-B'(k) upon a 
transmitted signal from an n-th transmitter using the 
channel value h mn (q) and the soft decision transmitted 
symbol b' n (k) as follows: 
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(2-D 



0 



// = ; 



! 0 



r a 



ii (?) 



L hMHq) 



(0) 

hMN(q) J 



(0-1 ) 



5'(A)=[^' T (^+Q-1) . . . fc' T (£) . . . fc /r 0H2+l)J r 
b'(k+q)=[b' 1 (k+q)b' 2 (k+q) . . . b' N (k+q)Y 
q=Q-l . . . -Q+l q*0 

b'(k)=\p\(k) . . . -ftb'jfi) . . . fe' N (*)] T ?=o 

where b'(k) has an element f(b' n (k)) at an n-th position, f( 
) is a function of a variable b' n (k) which satisfies that of 
f(0)=0 and is d{f(b' n (k))}/d{ b' n (k)}^0 and [] T is a 
transposed matrix; 

subtracting the interference component HB'(k) from a 
received matrix y(k) to define a difference matrix y'(k) 
where 

><£)=[r T (£+G-l)r T (A:+Q-2) . . . r^{k)f 



r(k)=[r ± (k)r 2 (k) 



defining a covariance matrix for noise components in the 
received matrix y(k) as cr 2 I which is determined from 
a variance a 2 of a Gaussian distribution and unit matrix 
I, and filtering the difference matrix y'(k) with an 
adaptive filter coefficient w n (k) which is determined by 



h = 



ffl,{Q-l)N+n 
H2,{Q-l)N+n 

HM-Q,{Q-l)N+n 



to provide a log-likelihood ratio for the received signal 
which corresponds to the transmitted signal from the 
n-th transmitter and from which interferences are elimi- 
nated. 

11. A turbo -reception method of receiving a transmitted 
signal from N transmitters where N is an integer equal to or 
greater than 2 comprising the steps of 



determining a channel value which is a transmission 
characteristic of a received signal from M received 
signals (where M is an integer equal to or greater than 
f) and a known signal;. 

estimating a soft decision transmitted symbol each from N 
a priori information; 

dividing N transmitted signals each into L (L^N) groups 
of transmitted signals including one or a plurality of 
transmitted signals, using a soft decision transmitted 
symbol and a channel matrix comprising channel val- 
ues to determine L equalized signals from which inter- 
ferences from other group of transmitting signals are 
eliminated and post -equalization channel information 
each corresponding to the transmission characteristic of 
the equalized signal for each group of transmitted 
signals; 



for each combination of L equalized signals and associ- 
ated channel information, treating the group of equal- 
ized signals as a received signal having a channel value 
defined by the channel information and where there is 
a plurality of transmitted signals which constitute the 
group of equalized signals, dividing such transmitted 
signal into a plurality of sub-group of one or a plurality 
of transmitted signals, using the soft decision transmit- 
ted signal to determine an equalized signal for the 
subgroup of transmitted signals from which interfer- 
ences from the other groups of transmitted signals are 
eliminated and associated post-equalization channel 
information, and where the group comprises a single 
transmitted signal, using the resulting equalized signal, 
the channel information and the soft decision transmit- 
ted signal to eliminate an interference of the transmitted 
signal itself due to multipaths; 

repeating the steps of dividing, eliminating interferences 
and generating post-equalization channel information 
until transmitted signals which constitute every equal- 
ized signal becomes a single one, thus finally deter- 
mining a equalized signal for each transmitted signal 
from which an interference caused by multipaths of its 
own is eliminated or determining, for the combination 
of the equalized signal and associated channel infor- 
mation, an equalized signal from which an interference 
between different transmitted signals and an intersym- 
bol interference by the transmitted signal itself is 
eliminated for each transmitted signal which consti- 
tutes the equalized signal. 

12. A turbo-reception method according to claim 11, in 
which for each group of transmitted signals, using the soft 
decision transmitted symbol and the channel value to gen- 
erate a replica of interferences caused by other groups of 
transmitted signals, subtracting the replica of interferences 
from the received signal to define a difference signal, 
determining, for each difference signal, a filter characteristic 
which is used to eliminate residual interference and associ- 
ated post-equalization channel information from the channel 
value and the soft decision transmitted symbol, filtering a 
corresponding difference signal with the residual interfer- 
ence eliminating filter characteristic to derive the equalized 
signal. 

13. A turbo-reception method according to claim 12, 
comprising the steps of 



(k), 



determining from the received signals r 1 (k), . 
a received matrix as defined below 

y(k)=[^(k+Q-iy^(k+Q-2) . . . r T (k)Y 
r(k)=[ ri (ky 2 (k) . . . r M (£)] T 
where [] T represents a transposed matrix; 

defining the transmission characteristic as a channel 
matrix H as defined below 



H = 



(0) 



0 



H iQ-i) 



H(0) 



0 



H (Q-i) 



hn(q)=[hm(<2) ■ ■ ■ KiJ({)Y 
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where m= 1 , . . . , M; n=l , . . . , N; and q=(), . . . , Q-l and 
Q is a number of muti -paths, and h mn (q)is a channel 
value of a path q from an n-th transmitter which is 
contained in the received signal r m ; 

denoting the soft decision transmitted symbol by b' n (k), 
calculating a replica HB'(k) upon one of groups of 
transmitted signals, formed by a first to a U-th trans- 
mitted signal where U is an integer satisfying N>U^1 
caused by other groups of transmitted signals as fol- 
lows; 

BXkyib'^k+Q-l) . . . b'^(k) . . . b' T (k-Q+l)Y 

b'(k+q)=[b , 1 (k+q)b' 2 (k+q) . . . b'Jk+q) . . . 2>' N (*+?)] 
T .?=0-l, ... 1 



b'(k+q)=[Q 



0 b n+1 '(k+q) . . . b'^{k+q)Y:q=0, 



where b'(k+q) have zero elements which are equal to U in 
number; 

subtracting the interference replica H-B'(k) from the 
received matrix y(k) to define the difference matrix y' 

(k). 

14. A turbo -reception method according to claim 11 in 
which when a further elimination of interferences is desired 
for the equalized signal and associated channnel informa- 
tion, the number of multipaths which are used during the 
elimination of interferences from the equalized signal is 
reduced. 

15. A turbo -reception method according to claim 14, 
further comprising the steps of determining from the 
received signals r 1 (k), . . . , r M (k), a received matrix defined 
as follows; 

y(k)=[r T (k+Q-l)r T (k+Q-2) . . . r T (A)J r 
r(k)=[rjk)r 2 (k) . . . r M (*)] T 

where [] T represents a transposed matrix; 

defining the transmission characteristic in the form of a 
channel matrix H, defined as follows; 



H = 



Hi 



(0) 



H, 



(0-D 



0 



(0) 



0 



Hio- 



(0-D 



H{q)=\h ± {q) . . . h N (qy] 

fcn(?)=l>ln(<?) • • - ^Mn(?)F 



where m=l, . . . , M; n=l, . . . N; and q=(), . . . , Q-l and 
Q is a number of mutipaths, and h mn (q)is a channel 
value of a path q from an n-th transmitter which is 
contained in the received signal r r 



denoting the soft decision transmitted symbol by b' n (k), 
one of the groups of transmitted signals comprising a 
first to U-th transmitted signal where U is an integer 
satisfying the inequality N>Ui^l, an equalized signal 
being formed for this group of transmitted signals so 
that interferences therein are eliminated by considering 
a number of multipaths equal to Q' where Q'<Q, 
calculating a replica of interference upon this group 
from other groups of transmitted in the form of HB'(k) 
as follows; 



B'(k)=[b' T (k+Q-l) . . . b'^k) . . . b^k-Q+l)? 

b'(k+q)=\b'Jk+q)b'Jk+q) . . . b'Jk+q) . . . b' N (k+q)] 
T :q=Q-l, ... 1 

b'(k+q)=[0 ... 0 b' u+1 (k+q) . . . b^(k+q)Y:q=0, . . . 

where b'(k+q) have 0 elements which are equal to U in 
number, and 

b'(k+q)=[b' 1 (k+q) . . . b'Jk+q) . . . b'^(k+q)Y :q=Q', . . 

■ -G+l 

subtracting the interference replica H-B'(k) from the 
received matrix y(k) to define the difference matrix y' g 

(k>. 

16. A time-reception method according to one of claims 1, 
3, 9, 10 and 11 in which during a second and a subsequent 
iteration of the turbo-reception processing, both the known 
singnal and the transmittes symbol hard decision output 
which is obtained during a previous alteration are used as 
reference signals, and the reference signals and the received 
signal are used to calculate the channel matrix. 

17. A turbo-reception method according to claim 16 in 
which one of the transmitted symbol hard decision outputs 
which are obtained during a previous iteration and which has 
a certainty in excess of a given value is also used as a 
reference signal to be used in the calculation of the channel 
matrix. 

18. A turbo -reception method according to one of claims 
1, 3, 9, 10 and 11 in which N a priori information X 2 [b n (k)] 
are derived from N decoders which correspond to the N 
transmitters, a log-likelihood ratio which is obtained for the 
received signal which corresponds to the n-th transmitted 
signal and through which interferences are eliminated is fed 
to a corresponding one of the decoders. 

19. A turbo -reception method according to one of claims 

1, 3, 9, 10 and 11 in which the N transmitted signals are 
transmitted from N transmitters which transmit a single 
information series in the form of N parallel series, the N a 
priori information information X 2 [b n (k)] being a result of 
series-parallel converstion of a priori information X 2 [b(j)] 
from one of the decoders, N log-likelihood ratios which are 
received signals corresponding to the N transmitted signals 
and from which interferences are eliminated being subject to 
a parallel-series convertion before being fed to the decoders. 

20. A turbo-receiver for receiving a signal from N trans- 
mitters where N represents an integer equal to or greater than 

2, comprising 

a received signal generator for forming M received sig- 
nals r m (where M is an integer equal to or greater than 
1) where m=l, . . . , M; 

a channel estimator to which each received signal r m and 
a reference signal in the form of a known signal are 
input for calculating a channel value h mn (q) and a 
channel matrix H, as defined below 



11 = 



H, 



(0) 



Hin- 



(0-D 



0 



H 



(0) 



0 



Hin- 



(0-D 



_h-Ml{q) ■■■ hMN(q) 
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-continued 



n = 1, 



input to calculate a channel value h mn (q) and a channel 
matrix H, as defined below. 



a received matrix generator which receives each received 
signal r m for generating a received matrix, as defined 
below 

r(k)=[r ± (k)r 2 (k) . . . r M (*)] T 

where k represents a discrete time, Q a number of 
multipaths of each transmitted wave, q=0, . . . , Q— 1, 
and [] T representing a transposed matrix; 

a soft decision symbol generator receiving N a priori 
information for generating a soft decision transmitted 
symbol b' n (k); 

a replica matrix generator to which respective soft deci- 
sion transmitted symbols b-^k) to b N '(k) are input to 
generate an interference replica matrix B'(k) with 
respect to an n-th transmitted signal as follows; 

B'(*)=[£' T (fc+!2-l) . . . b'^ik) . . . b' T (k-Q+l)Y 
bXk+q)=\b' ± {k+q)b'^k+q) . . . b' N (k+q)Y 
q=Q-l, . . . ,-Q+l q*0 
b'(k)=[b\(k) . . .0 . . . b' N (k)Y q=0 

where b'(k) has a zero element at n-th position; 

a filter processor to which the channel matrix H and the 
interference replica matrix B'(k) are input for calculat- 
ing and delivering an interference component H-B'(k) 
with respect to a received signal corresponding to the 
n-th transmitted signal; 

a difference calculator to which the interference compo- 
nent HB'(k) and the received matrix y(k) are input for 
delivering a difference matrix y , (k)=y(k)-H-B , (k); 

a filter coefficient estimator to which the channel matrix 
H or a reference signal is input for determining an 
adaptive filter coefficient w n (k)to be applied to a 
received signal which corresponds to a transmitted 
signal from the n-th transmitter in order to eliminate 
residual interference components which remain in the 
reference matrix y'(k); 

and an adaptive filter to which the difference matrix y'(k) 
and the adaptive filter coefficient w n (k) are input to 
filter y'(k) to provide a log-likelihood ratio as an 
received signal which corresponds to the transmitted 
signal from the n-th transmitter and from which inter- 
ference is eliminated and which is then fed to the n-th 
decoder. 

21. A turbo-receiver for receiving a signal from N trans- 
mitters where N is an integer equal to or greater than 2, 
comprising 

a received signal generator for forming M received sig- 
nals r m (where M is an integer equal to or greater than 
1) where m=l, . . . , M; 

N decoders; 

a channel estimator to which each received signal r m and 
a reference signal in the form of a known signal are 



// = 



(0) 



(2-1) 



0 



(0) 



0 



n = 1, .. . N 



n=l, . . . , N 

a received matrix generator to which each received signal 
r m is input to generate a received matrix, as defined 
below 

y(k)=[r T (k+Q-iy T (k+Q-2) . . . r T (k)Y 
r(k)=[ ri (k)r 2 (k) . . . r M (£)Y 

where k represents a discrete time, Q a number of 
multipaths of each transmitted wave, q=0, . . . Q-l and 
[] T a transposed matrix; 

a soft decision symbol generator to which N a priori 
information are input to generate soft decision trans- 
mitted symbols b' n (k) (where n=l, . . . , N); 

a replica matrix generator to which the soft decision 
transmitted symbols b\ (k) to b' n (k) are input to 
generate an interference replica matrix B'(k) for a 
transmitted signal from an n-th transmitter as indicated 
below 

B\k)=[b' T (k+Q-l) . . . b'^ik) . . . £' T (/t-<2+l)F 
bXk+q)=[b\(k+q)b f 2 (k+q) . . . b'^(k+q)Y 
q=Q-l, . . . ,-Q+l q*0 

b'(k)=[b\(k) . . . -fib'Jky) . . . fc' N (*)] T q=0 

where b'(k) has an element f(b' n (k)) at an n-th position, 
and f( ) is a function of a variable b' n (k) which satisfies 
that f (0)=0 and 

d{f(b* n (k))}/{ b' n (k)}^0; 

a filter processor to which the channel matrix H and the 
interference replica matrix B'(k) are input to calculate 
and deliver an interference component H-B'(k) for a 
received signal which corresponds to the transmitted 
signal from the n-th transmitter; 

a difference calculator to which the interference compo- 
nent H-B'(k) and the received matrix y(k) are input to 
deliver a difference matrix y'(k)=y(k)-H-B'(k); 

a filter coefficient estimator to which the channel matrix 
H or a reference signal is input to determine an adaptive 
filter coefficient w n (k) for a received signal which 
corresponds to a transmitted signal from the n-th trans- 
mitter in order to eliminate residual interferences in the 
difference matrix y'(k); 

and an adaptive filter to which the difference matrix y'(k) 
and the adaptive filter coefficient w n (k) are input to 
filter y'(k) to provide a log-likelihood ratio as a received 
signal which corresponds to the transmitted signal from 
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the n-th transmitter and from which interferences are 
eliminated and which is then fed to the n-th decoder. 

22. A turbo -receiver for receiving a transmitted signal 
from N transmitters where N is an integer equal to or greater 
than 2, comprising 

a received signal generator for forming M received sig- 
nals (where M is an integer equal to or greater than 1); 

a channel estimator to which the N received signals and 
a reference signals in the form of a known signal are 
input to estimate a channel value representing a trans- 
mission characteristic; 

a prestage equalizer to which the M received signals, the 
channel value and N a priori information are input for 
delivering, for each transmitted signal from one or 
more transmitters, a plurality of sets of equalized signal 
from which interferences by transmitted signals from 
other transmitters are eliminated and an associated 
post-equalization channel information; 

and a plurality of poststage equalizers to which the set of 
the equalized signal and associated channel informa- 
tion from the prestage equalizer, and a priori informa- 
tion which corresponds to the transmitted signal which 
constitutes the equalized signal are input to deliver a 
log-likelihood ratio by eliminating from the equalized 
signal an intersymbol interference from each transmit- 
ted signal which constitutes the equalized signal due to 
the multiple paths and a mutual interferences between 
each transmitted signal and other transmitted signals 
which also constituted the equalized signal. 

23. A turbo -receiver for receiving a transmitted signal 
from N transmitters where N is an integer equal to or greater 
than 2, comprising 

a received signal generator for forming M received sig- 
nals (where M is an integer equal to or greater than 1); 

a channel estimator to which the N received signals and 
a reference signal in the form of a known signal are 
input to estimate a channel value representing a trans- 
mission characteristic; 

a prestage equalizer to which the M received signals, the 
channel value and N a priori information are input to 
deliver, for each transmitted signal from one or more 
transmitters, a plurality of sets of an equalized signal 
from which interferences by transmitted signals from 
other transmitters are eliminated and an associated 
post-equalization channel information; 

and a plurality of poststage equalizers to which the sets of 
the equalized signal and associated channel informa- 
tion from the prestage equalizer, and a priori informa- 
tion which corresponds to a plurality of transmitted 
signals which constitute the equalized signal are input 
to deliver, for each or for a sub-group of transmitted 
signals among the plurality of transmitted signals 
which constituted the equalized signals, a plurality of 
sets of an equalized signal from which interference by 
other transmitted signals which constituted that the 
equalized signals are eliminated and an associated 
post-equalization channel information. 

24. A turbo- receiver according to one of claims 20 to 23, 
further comprising a previous symbol memory in which a 
hard decision transmitted symbol from the decoder is stored 



to be updated, and means for reading the hard decision 
transmitted symbol from the previous symbol memory to 
feed it to the channel estimator as the reference signal during 
a second and a subsequent iteration of the turbo-reception 
processing. 

25. A turbo -receiver according to claim 24 further com- 
prising a comparator for comparing a soft decision trans- 
mitted symbol which is input thereto against a threshold 
value, and a selector which is controlled by an output from 
the comparator so that one of the hard decision transmitted 
symbols, for which a corresponding soft decision transmit- 
ted symbol has a value in excess of the threshold value, is 
stored in the previous symbol memory. 

26. A turbo-receiver according to one of claims 20 to 23, 
further comprising N decoders, to which the N log-likeli- 
hood ratios are delivered, the N decoders providing the N a 
priori information at their outputs. 

27. A turbo-receiver according to one of claims 20 to 23 
in which the N transmitted signals are transmitted signals as 
the N transmitters transmit a single information series as N 
parallel series, further comprising a parallel- series converter 
for converting the delivered N log-likelihood ratios into 
serial series, a decoder to which the log-likelihood ratio in 
the serial series is input, and a series-parallel converter for 
converting a priori information from the decoder into N 
parallel series to provide the N a priori information. 

28. A turbo-reception method in which a channel value 
representing a transmission path characteristic of a received 
signal is estimated from the received signal and a known 
signal serving as a reference signal, the received signal is 
processed in accordance with the estimated channel value, 
the processed signal is decoded, and the processing which 
uses the estimated channel value and the decoding are 
iterated upon the same received signal; farther comprising 

determining the certainty of a decoded hard decision 
information signal on the basis of a value of a corre- 
sponding soft decision information symbol, and using 
a hard decision information symbol having a certainty 
which is in excess of a given value in the channel 
estimation of a next iteration as a reference signal. 

29. A reception method according to claim 28, further 
comprising the step of calculating a 2 I as a covariance 
matrix of noise component in a received matrix y(k) where 
a 2 represents a variance of Gaussian distribution and I a unit 
matrix. 

30. A reception method according to claim 28 further 
comprising using an estimated channel matrix H and a 
received signal matrix y(k) to calculate a covariance matrix 
U of noise components within the received signal matrix 
y(k) during each iteration as follows; 

£(£)=|> T (A:+e-l) . . . b^k) . . . b^k-Q+^Y 

b(k+q)=\bjk+q) ■ ■ ■ & N (*+<7)J T (q=-Q+i ■ ■ ■ e-i) 

b 1 (k+q) to b N (k+q) being the reference signals compris- 
ing the known signal and the hard decision information 
symbols having a certainty which is in excess of a given 
value, Tr being the length of the reference signal. 

31. A reception method according to claim 28 in which the 
iteration of the processing which uses the estimated channel 
value and decoding processing comprises an iteration of 
determining a linear equalization filter in accordance with 
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the estimated channel value, processing the received signal 
by the linear equalization filter and decoding the processed 
signal. 

32. A reception method according to claim 28 in which the 
iteration of the processing which uses the estimated channel 
value and decoding processing comprises an iteration of 
performing a rake synthesis which compensates for a phase 
rotation which each symbol has experienced on a transmis- 
sion path in a rake synthesis processor in accordance with 
the estimated channel value, and decoding the signal which 
is produced by the rake synthesis in a turbo -de coder. 

33. A reception method according to claim 28 in which the 
iteration of the processing which uses the estimated channel 
value and decoding processing comprises an iteration of 
setting up weights which determine an antenna directivity 
response of an adaptive array antenna receiver in accordance 
with the estimated channel value, and decoding an output 
from the adaptive array antenna receiver in a turbo -decoder. 

34. A reception method according to claim 33 in which a 
rake synthesis which compensates for a phase rotation which 
each symbol has experienced on a transmission path in 
accordance with the estimated channel value is made for an 
output from the adaptive array antenna receiver in a rake 



synthesis processor, and a resulting signal from the rake 
synthesis is fed to the turbo-decoder. 

35. A receiver in which a channel value representing a 
transmission path characteristic of a received signal is 
estimated from the received signal and a known signal 
serving as a reference signal, the received signal is processed 
by using the estimated channel value, the processed signal is 
decoded, and the processing which uses the estimated chan- 
nel value and the decoded processing are iterated upon the 
same received signal; further comprising 

means for determining whether or not the certainty of a 
decoded hard decision information symbol is in the 
excess of a given value by seeing whether or not a 
corresponding soft decision information symbol has a 
value which is in excess of a threshold value; 

and a previous symbol memory having its stored content 
updated by a hard decision information symbol which 
has been determined as likely to be certain, the stored 
content of the previous symbol memory being used as 
a reference signal during the channel estimation of a 
next ilteration. 



